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Zusammenfassung 


Die Formulierung und Implementierung neuartiger Materialmodelle 
auf der Mikroskala polykristalliner Materialien, z.B., kubisch fláchen- 
zentrierte Metalle, gibt ein tiefgründiges Verständnis inner- und inter- 
kristalliner Effekte. Dies erlaubt eine genauere Vorhersage physikali- 
scher Effekte, welche das Verhalten des Materials auf der makrosko- 
pischen Skala beeinflussen. Im Kontext der erweiterten Kontinuums- 
mechanik wird eine Plastizitátstheorie basierend auf physikalischen 
Effekten von Versetzungen hergeleitet. Die erforderlichen Annahmen 
werden konstruiert um eine Theorie äquivalent zu einem Modellansatz 
basierend auf dem Prinzip virtueller Kráfte zu erhalten. Es wird eine 
orientierungsabhängige Fließbedingung der Korngrenzen entwickelt. 
Mittels einer analytischen Lösung eines dreiphasigen periodischen 
Laminates wird die Aufstauung der Versetzungen an Korngrenzen 
und der Transport über die Korngrenzen hinweg untersucht. Dies 
ist von besonderem Interesse bei der Untersuchung der plastischen 
Verformung von Metallen. Der sich ergebende Zusammenhang der 
Festigkeit und der Korngröße ist konsistent mit der Materialtheorie 
der Versetzungen. Im Gegensatz zum üblich verwendeten Ansatz der 
Fließbedingung der Korngrenzen besteht der neu entwickelte Ansatz 
die Konsistenzprüfung der Einkristalllösung für den Grenzfall aufein- 
andertreffender Gleitsysteme. Für missorientierte Körner ergibt sich ein 
diskontinuierlicher Verlauf der plastischen Abgleitung, was eine beson- 
dere Behandlung in der Anwendung mittels Finiter Elemente erfordert. 
Die Fließbedingung der Korngrenzen wird an scharfen Grenzschich- 
ten mittels diskontinuierlicher Ansatzfunktionen ausgewertet. Hierfür 
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wird ein physikalisch vereinfachtes, numerisch effizientes Modell unter 
Annahme großer Deformationen eingeführt. Um die weitreichenden 
Interaktionen der Versetzungen zu erfassen, enthält das Modell eine 
Gradientenspannung, welche auf den Gradienten der akkumulierten 
plastischen Abgleitung beschränkt ist. Es wird ein Ansatz formuliert 
um den rechnerischen Nutzen akkumulierter Feldvariablen in einem 
physikalisch motivierten Framework basierend auf dem Transport 
von Versetzungen aufrecht zu erhalten. Mittels anteiliger Verteilung 
des planaren Gradienten der akkumulierten plastischen Abgleitung 
auf die einzelnen Gleitsysteme wird ein numerisch effizientes Frame- 
work für die mechanismusbasierte Modellierung der Kristallplastizität 
konstruiert. Der Einfluss der Modellparameter auf die mechanische 
Antwort einer Laminatstruktur und einer Faser-Matrix Struktur wird in 
dreidimensionalen Finiten Elemente Simulationen untersucht. Die sich 
ergebenden Verteilungen der plastischen Abgleitung und der Verset- 
zungsdichte sind konsistent mit experimentellen Untersuchungen und 
zeigen die Notwendigkeit auf, Energieterme und Fließregeln basierend 
auf Versetzungstheorien für die Modellierung der Interaktionen von 
Versetzungen und Korngrenzen zu nutzen. 
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Summary 


For polycrystalline materials, e.g., face-centered cubic metals, the for- 
mulation and implementation of novel material models on the mi- 
croscale gives a deeper understanding of inner- and inter-crystalline 
effects. This allows more precise prediction of physical effects which 
influence the material on the macroscopic scale. A physically-based dis- 
location theory of plasticity is derived within an extended continuum 
mechanical context. The required assumptions are constructed in order 
to obtain an equivalence to a model approach based on the principle 
of virtual power. An orientation dependent grain boundary flow rule 
is introduced. With an analytical solution of a three-phase periodic 
laminate the dislocation pile-up at grain boundaries and transmission 
through the grain boundaries is investigated, which is of special inter- 
est during the plastic deformation of metal materials. The observed 
dependence of the material strength on the grain size is consistent 
with underlying dislocation theories. In contrary to the common 
approach of grain boundary modeling the introduced grain boundary 
flow rule passes the single crystal consistency check for the limit case 
of coinciding slip systems. For misaligned grains a discontinuous 
slip distribution is obtained, which requires special treatment in the 
finite element application. The grain boundary flow rule is evaluated 
at sharp interfaces using discontinuous trial functions. Hereby, a 
physically simplified, numerically efficient approach is introduced in 
the finite strain setting. The model is restricted to one gradient-stress, 
associated with the gradient of the accumulated plastic slip, in order 
to account for long range dislocation interactions. An approach is 
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Summary 


formulated in order to regain the computational benefit of accumulated 
plastic field variables within a physically motivated framework based 
on dislocation transport. By prorating the in-plane gradient of the 
accumulated plastic slip on the individual slip systems a numerically 
efficient framework for mechanism based modeling of crystal plasticity 
is constructed. The influence of model parameters on the mechanical 
response of a laminate grain structure and a fiber-matrix composite 
material is investigated in three-dimensional finite element simulations. 
The evolving distributions of plastic slip and dislocation density are 
consistent to experimental observation and show the importance of 
energy expressions and flow rules based on underlying dislocation 
theories for the modeling of dislocation-grain boundary interactions. 
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Introduction 


11 Motivation and research objectives 


In this work the dislocation-induced plastic deformation of polycrys- 
talline face-centered cubic (FCC) metals is considered. During the plas- 
tic deformation dislocations arise inside the bulk material and are trans- 
ported on slip systems through the inner crystalline structure. The 
transport of dislocations, however, is hindered by the presence of grain 
boundaries (GBs). The pile-up of dislocations at and their transition 
across GBs has a significant impact on the overall material behavior. 
Therefore, the modeling of the dislocation-GB interactions is of utmost 
interest. Discrete dislocation dynamics or molecular dynamics simula- 
tions can be performed to investigate the movement of dislocations and 
interactions of dislocations with obstacles or with each other. For the 
application of polycrystalline structures, however, these methods are 
not applicable due to the high computational cost. In a continuum me- 
chanical approach plastic slips, which result from the accumulation of 
several gliding dislocations, are considered. In crystal plasticity frame- 
works the mechanical behavior of the material is modeled with the 
construction of constitutive laws for the evolution of plastic slips based 
on underlying dislocation theories. 

Classical continuum theories, however, fail to reproduce dependency 
of the material strength on the grain size, which is observed experimen- 
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tally for polycrystalline materials. A phenomenological approach is 
given by the introduction of a grain size dependency of the slip system 
yield strength, i.e., the stress which is required to act on a slip system 
for the activation of plastic slip evolution. As concluded by Zhu et al. 
(2008), however, each of these models explains size effects in differ- 
ent strain regions, under different initial conditions and for different 
deformation mechanisms. In order to model a physically based size- 
dependent behavior of polycrystalline metallic materials, the plastic 
slips are introduced as additional degrees of freedom (DOF). A close 
connection exists between the gradient of plastic slip and the density of 
geometrically necessary dislocations (GND), cf. Ashby (1970). Since the 
early work of Aifantis (1987) many gradient plasticity frameworks have 
been developed. The formulation of models which are optimal with re- 
gard to the dispute between physical correctness and implementability 
is an active field of current research. Besides applications in finite ele- 
ment (FE) simulations analytical solutions can be constructed to clarify 
the influence of energy expressions and dislocation-GB interactions. 

Within the thesis, several originalities were developed, whereas the 

main highlights are: 

* The construction of a gradient crystal plasticity framework based on 
an extended balance equation. The comparison to a model approach 
based on the principle of virtual power is performed. The required as- 
sumptions are formulated in order to obtain an equivalence between 
both approaches. 

* The formulation of a GB flow condition based on the geometric dis- 
location tensor. The flow condition couples plastic slips at GBs by a 
misorientation function of adjacent slip system. The flow condition 
describes the orientation dependent resistance of a GB. This allows 
to model the pile-up of dislocations and the activation of GB slip, i.e., 
the transmission of dislocations across the GB. 
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* The construction of an analytical solution for the distribution of plas- 
tic slip and the displacement field of a periodic three-phase laminate 
structure. The analytical solution is used for a comparison of the 
GB flow condition approach to a common approach based on the 
assumption of independent slips at the GB. The analytical solution 
gives an insight into the GB behavior for misaligned grains, the con- 
sistency check for the limiting case of a vanishing misorientation and 
the investigation of size effects. 

* The application of shape function enrichment in the FE implementa- 
tion of gradient crystal plasticity frameworks. This method allows 
for discontinuous solution fields at interfaces and is well-suited for 
the modeling of dislocation transmission across GBs. 

* The extension of a numerically efficient gradient crystal plasticity 
framework based on the accumulated plastic slip to finite strains and 
discontinuous solution fields. The model accounts for the pile-up of 
dislocations at GBs and the transmission of dislocations across GBs 
by a GB energy dependent on the jump and mean value of plastic 
slip. The FE implementation provides insight into the GB behavior 
of a periodic laminate structure and the investigation of material 
parameter influences. 

* The development and implementation of a finite strain gradient crys- 
tal plasticity framework that accounts for a semi-phenomenological 
dislocation transport mechanism. A physically motivated disloca- 
tion transport equation is constructed while preserving the compu- 
tational benefit of accumulated field variables. The model is based 
on a dislocation energy which is obtained by a statistical theory of 
dislocations. This allows for a semi-phenomenological treatment of 
dislocations within the context of gradient crystal plasticity. 
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1.2 Outline of the thesis 


In Chapter 2 an introduction of dislocations and dislocation-based plas- 
ticity is given. Fundamental informations on the origin of dislocations 
and strengthening mechanisms associated with dislocation multiplica- 
tion and dislocation-GB interactions are briefly presented. 

In Chapter 3, fundamentals of classical continuum mechanics such as 
kinematics and basic assumptions of small strain and finite strain frame- 
works are presented. Based on mechanical and thermodynamical bal- 
ance relations for a material divided by a singular surface the Clausius- 
Duhem inequality is obtained for regular and singular points. 

This provides the basics for Chapter 4, in which a gradient crystal plas- 
ticity theory is derived by an extension of balance equations. Potential 
relations and boundary conditions for both, bulk material and GBs, are 
obtained. A comparison to a model approach based on the principle of 
virtual power is made. Simplifying assumptions and their implication 
on the Clausius-Duhem inequality and the bulk material flow rule are 
presented in order to obtain consistent frameworks. 

Chapter 5 introduces a GB flow condition approach. The GB energy is 
introduced as a function of the geometric dislocation tensor in order to 
account for the misorientation of adjacent slip systems and dislocation 
transmission across GBs. An analytical solution for the distribution of 
plastic slip and the displacement field of a periodic laminate is derived. 
For simplicity small strains and single slip is considered. The GB flow 
condition approach is compared with a commonly used approach for 
a varying grain misorientation angle. Consistency checks are made for 
the limiting case of vanishing grain misorientation. In addition, the 
influence of the grain size is investigated. 

In Chapter 6 a finite strain gradient crystal plasticity theory is presented, 
which includes a FE application based on a micromorphic approach. 
In order to reduce the computational cost the model is restricted to 
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one gradient-stress associated with the gradient of the accumulated 
plastic slip. The Clausius-Duhem inequality and flow rules for the bulk 
material and GBs are derived. Details on the FE implementation of a 
discontinuous plastic slip field are presented. The influence of grain 
misorientation, grain size and material parameters are discussed for a 
periodic laminate microstructure. 

Chapter 7 introduces a finite strain gradient crystal plasticity frame- 
work that accounts for a semi-phenomenological dislocation transport 
mechanism. Besides the accumulated plastic slip an accumulated 
density of dislocations is introduced as an additional degree of freedom. 
Flow rules for the bulk material and the GBs are derived. Three- 
dimensional FE simulations of a periodic laminate microstructure and 
a fiber-matrix composite material are performed. 

A summary in Chapter 8, featuring an outlook based on remaining 
open questions from the work at hand, concludes the thesis. 


13 Frequently used acronyms, symbols, and 


operators 
Acronyms 
BCC Body-centered cubic 
DOF Degree of freedom 
FCC Face-centered cubic 
FE Finite element 
GB Grain boundary 
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GND Geometrically necessary dislocations 


SSD Statisticall stored dislocations 


Latin letters 


a,b, A,B,D,... Scalar quantities 

U,V,Ww,... First-order tensors 

A,B,C,... Second-order tensors 

A,B,C,... Fourth-order tensors 

A, Micro-inertia force 

b Absolute value of the Burgers vector 

by Micro-body force 

b Body force 

B Finite element shape function derivative 

CRC Intra-grain interaction moduli 

GC Inter-grain interaction moduli 

C Fourth-order elastic stiffness tensor 

da Slip direction of slip system a 

D Dissipation 

Died Reduced dissipation 

e Specific internal energy 

€1, €2, €3 Unit vector in x-,y- and z-direction of a Cartesian 
coordinate system 

E, Green strain tensor 


DAFN 
"tj 
"S 


= 


32 T 


Deformation gradient 
Elastic and plastic part of F 
Geometric dislocation tensor 
Penalty modulus 
Displacement gradient 
Elastic and plastic part of H 
Second-order identity tensor 
Fourth-order identity tensor 
Mean curvature 

Isotropic energy parameter 
Correlation energy parameter 


Defect energy parameter 


1.3 Latin letters 


Grain boundary energy parameter associated 


with the jump of plastic slip 


Grain boundary energy parameter associated 


with the mean value of plastic slip 


Dislocation self energy parameter 


Velocity gradient 
Distortion rate tensor 


Schmid tensor of slip system a 


Symmetric Schmid tensor of slip system a 


Normal vector 


Slip plane normal vector of slip system a 


Finite element shape function 


Slip rate sensitivity parameter 
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Generalized relative stress 
Projector of slip system a 
Projector of interface S 
External power 

Internal power 

Heat flux vector 

Orthogonal rotation matrix with angle ù 
Inter-grain rotation matrix 
Residuals 

Rigid body rotation tensor 
Piola-Kirchhoff stress 
Dislocation direction line 
Second Piola-Kirchhoff stress 
Time 

Stress vector 

Surface traction 
Displacement field 

Rigid stretch tensor 

Velocity field 


Position vector of a material point in the current 


configuration 


Position vector of a material point in the refer- 


ence configuration 


Greek letters 


Ba 


1.3 Greek letters 


Slip system index 


Derivative of hardening energy by plastic slip of 
slip system a 


Derivative of hardening energy by dislocation 
density 


Plastic slip of slip system a 
Accumulated plastic slip of slip system a 
Accumulated plastic slip 
Micromorphic accumulated plastic slip 
Infinitesimal strain tensor 

Elastic and plastic part of € 

Continuous plastic slip of slip system o 
Specific entropy 

Temperature 

Initial hardening modulus 

Saturation hardening modulus 
Orientation angle of single slip system 
Grain boundary strength 

Plastic slip parameter of slip system a 
Reference shear rate 

Thermal conductivity 

Shear modulus 


Dislocation velocity 
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Coarse grained dislocation velocity 

Grain boundary misorientation parameter 
Vector-valued microstress 

Microscopic surface traction 


Microscopic surface traction work conjugated to 
the jump and mean value of plastic slip 


Scalar-valued microstress 
Microtraction 

Density of dislocations 

Coarse grained density of dislocations 
Density of geometrically necessary dislocations 
Accumulated density of dislocations 
Initial density of dislocations 

Mass density 

Cauchy stress tensor 

Resolved shear stress 

Initial critical resolved shear stress 
Saturation stress 

Drag stress 

Yield function 

Specific free energy 

Specific correlation energy 

Specific elastic energy 


Specific defect energy 
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Specific hardening energy 
Specific dislocation self energy 
Specific numerical penalty energy 


Heat supply density 


Absolute value of a scalar quantity 

Norm of a vector or tensor 

Lagrangian description of a quantity 

Coarse graining of a quantity 

Quantity with respect to the interface S 
Jump of a quantity across an interface 

Mean value of a quantity across an interface 
Left-hand limit of a quantity at an interface 
Right-hand limit of a quantity at an interface 
Macaulay bracket: ramp function of a quantity 
Quantity in local coordinate system 

Material time derivative of a quantity 

Dot product of two tensors A, B 

Dyadic product of two tensors A, B 


Box product of two tensors A,B with 
(ADB)[C] = ACB 


Linear mapping of a second-order tensor 
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Linear mapping of a second-order tensor by a 


fourth-order tensor 

Inverse of a tensor A 

Transpose of a tensor A 

Cosine of a scalar quantity 

Determinant of a tensor 

Variation of a quantity 

Partial derivative of a quantity w.r.t., e.g., x 
Total derivative of a quantity w.r.t., e.g., x 
Eulerean divergence of a vector or tensor 
Surface divergence of a vector or tensor 
Lagrangian divergence of a vector or tensor 
Logarithm of a scalar quantity 

Eulerean gradient of a quantity 

Surface gradient of a quantity 

Lagrangian gradient of a quantity 
Lagrangian in-plane gradient of a quantity 
Sign of a scalar quantity 

Sine of a scalar quantity 

Symmetric part of a tensor 


Trace of a tensor 
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Introduction to dislocations 


2.1 Dislocation-based plasticity 


This chapter serves as an introduction to dislocations and dislocation- 
based plasticity. It gives a short overview of well-established theories 
and does not contain new findings. For a more detailed introduction to 
dislocations it is referred to Hull and Bacon (2011). 

An important characteristic of a dislocation is the Burgers vector. A 
Burgers vector can be constructed by creating a closed atom-to-atom 
loop around a dislocation. After reconstructing the same loop in a dis- 
location free crystal structure, the circuit shows a gap. The vector which 
results by completing the circuit, called Burgers vector b, determines the 
magnitude and orientation of the dislocation. With the Burgers vector 
and the dislocation direction line s it is possible to distinguish between 
an edge dislocation (s L b) and a screw dislocation (s || b). A graphical 
illustration of an edge dislocation in a crystal structure, marked by the 
symbol “L” is shown in Fig. 2.1. 

The discovery of dislocations by Orowan (1934), Polanyi (1934) and 
Taylor (1934) provides an explanation for the process of plastic defor- 
mation in crystals. Experimental investigations of Schmid and Boas 
(1935) show that due to dislocation motion the critical resolved shear 
stress of single crystals is found to be considerably smaller than the 
theoretical strength. The motion of dislocations can be classified into 
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Figure 2.1: Graphical illustration of an edge dislocation with Burgers vector b. 


conservative dislocation glide and non-conservative dislocation climb. 
Gliding of dislocations describes the motion on specific planes, which 
depend on the crystal structure of the material. Climb, in contrary, 
describes the motion out of these planes. In this work the focus is on 
dislocation glide. Glide of many dislocations results in slip, which is the 
most common manifestation of plastic deformation in crystalline solids, 
cf. Hull and Bacon (2011). The slip planes in FCC metals, in contrary 
to body-centered cubic (BCC) metals, are well defined as planes with 
the highest density of atoms. The direction of slip is the direction in 
which the atoms are most closely spaced. FCC crystals have four slip 
planes with three slip directions in each plane, and therefore twelve 
slip systems. Each slip system is characterized by the corresponding 
Schmid tensor Ma = da ® Nna, with the slip plane normal vector nq 
and the slip direction da, where a € {1,..., N} denotes the slip system 
index and N is the number of slip systems. The two unit vectors are 
orthogonal to each other, i.e., d: n, =0. In Appendix A.1 the slip 
system convention for FCC, which is used in this work, is given. The 
projection of the Cauchy stress tensor ø to the slip system a is intro- 
duced as the resolved shear stress 


Ta 7 0: M. (2.1) 
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2.2 Dislocation strengthening mechanisms 


As soon as the resultant shear stress reaches a critical value r. plastic 
slip occurs. The Orowan equation for plastic slip 


Ya = Pabva (2.2) 


gives a relation between the rate of plastic slip ^, and the average glide 
velocity of dislocations va in the slip system a. Hereby, b denotes the 
absolute value of the burgers vector, i.e., b = ||b|| and pa denotes the 
density of dislocations in the slip system «. 


22 Dislocation strengthening mechanisms 


Polycrystals consist of many grains with differently oriented slip sys- 
tems. Hence plastic flow first occurs in the grain with the slip system 
which can resolve the most shear stress. This leads to an inhomoge- 
neous plastic material behavior. Because the cohesion of the material 
still has to be guaranteed, the formation of GND, as illustrated in 
Fig. 2.2, occurs. Ensuing from the deformation of a polycrystal (a) 
where each grain behaves independently, gaps and overlaps would 
emerge (b). Due to dislocations these defects are compensated by shear 
(c) (top) or deformations normal to the GB (c) (bottom). By combining 
these processes the counteract of all possible failures is enabled (d). The 
interfaces between regions of different crystal orientations impede the 
motion of dislocations. The resistance of a GB against the transmission 
of dislocations is strongly influenced by the mismatch of adjacent slip 
systems. For an overview of basic dislocation mechanisms near GBs 
and transmission of dislocations across GBs it is referred Chapter 5 as 
well as to the review article of Bayerschen et al. (2016). Due to the 
pile-up of dislocations at GBs the grain size influences the mechanical 
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or DD 
: overlap : | "4 > 


(a) (b) (c) (d) 


Figure 2.2: Illustration of the origin of GND, cf. Ashby (1970). 


response of polycrystalline materials. The Hall-Patch relation 


Tc — TQ X 


(2.3) 


1 
Nx 
cf. Hall (1951), relates the difference between the yield stress 7. and 
the yield stress of an undeformed single crystal ro to the grain size d. 
This relation, however, is only a manifestation of the general size effect 
and can not be explained by a variety of theories found in literature, as 
discussed by Li et al. (2016). 

On the contrary to GND, statistically stored dislocations (SSD) are not 
required for a cohesion of the crystal and originate at various locations 
inside each grain. Essential origination processes are an unideal crystal 
growth and the Frank-Read-Mechanism illustrated in Fig. 2.3. A dis- 
location line lies in a slip plane, is pinned at both ends by two fixing 
points A and B, e.g., by crystal defects (a) and bows out due to an 
applied resolved shear stress 7. The radius of curvature R decreases 
as T increases until the minimum value of R is reached (b). As the 
line continues to expand, the dislocation becomes unstable and R in- 
creases (c). Consequently, this behavior leads to an annihilation of the 
opposing segments (d, e). The result is a large outer dislocation loop, 
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ey N 
L 
(d) (e) 


(a) (b) (c) 


Figure 2.3: Illustration of the Frank-Read mechanism, cf. Hull and Bacon (2011). 


which continues to expand, and a regenerated inner dislocation line, 
which repeats the process. The necessary stress for this Frank-Read 
mechanism is given with shear modulus G by 


ach 


LC (2.4) 


To = 
cf. Hull and Bacon (2011), where L denotes the spacing between the 
fixing points A and B and a z 0.5 denotes a material constant. A spe- 
cial case of the treatment above is the phenomenon of work hardening. 
Here, the pinning points are forest dislocations and the relevant size 
parameter is given by their average spacing L = 1/,/p, cf. Arzt (1998). 
The classical Taylor equation for work hardening is obtained 


Te = aGb/p. (2.5) 


Both types of dislocations, GND and SSD, influence the material behav- 
ior of polycrystals. At first an overview of classical continuum mechan- 
ics theory for the local continuum representation of dislocation-based 
plasticity is given. 
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Chapter 3 


Classical continuum mechanics 
theory 


3.1 Introduction and motivation 


The critical resolved shear stress evolves during the plastic deforma- 
tion process of the crystal lattice due to work-hardening. The work- 
hardening of single-crystals describes the strengthening of the mate- 
rial due to the production and annihilation of dislocations and is of- 
ten formulated by a constitutive relation between the critical resolved 
shear stress and scalar-valued hardening quantities, cf. , e.g., Bishop 
and Hill (1951); Hill (1966). With the formulation of constitutive laws 
for the evolution of the plastic slips, plastic deformation of polycrys- 
talline structures can be modeled. This classical crystal plasticity ap- 
proach is presented in the following. Hereby, mechanical and ther- 
modynamical balance relations are introduced in order to describe the 
motion of a material body in time and the evolution of the associated 
mechanical fields. 


3.2 Kinematics and basic assumptions 


The motion of a material body Y € R? is considered. As seen in Fig. 3.1, 
the position of a material point can be characterized by the vector X 
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fu) t0 


OVo 


€3 


Figure 3.1: Material body in the reference and the current placement. 


in reference description (time: t = 0), or by the vector x in spatial de- 
scription (time: t > 0). With the introduction of the deformation x the 
current position of a material point can be described as a function of its 
reference position upon neglecting phenomena like, e.g., cracks 


x= x(X,t). (3.1) 


It is assumed that the function x is invertible and twice continuously 
differentiable. The material time derivative of the motion represents 
the velocity v(X,t) = x(X,t). For a general field quantity (æ, t) the 
material time derivative is given as 

Ov(m.t) | Ov(a.t) 


w(a,t) = aE + — bz : v(a,t). (3.2) 


The displacement of a material point results as the difference of its 
current position and reference position 


u(X,t) = xX(X,t) - X. (3.3) 
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The displacement gradient reads 


H = Grad (u(X,t)) = u : (3.4) 
In the following, dependencies of quantities are dropped due to better 


readability. 
Small deformation setting 


In the context of Chapter 4 and 5 only small deformations are consid- 
ered. Frameworks extended to finite deformations are presented in 
Chapter 6 and 7. For a geometrically linear theory the norm of the 
displacement gradient has to fulfil the condition 


[H| «1, (3.5) 


and occurring rotations have to be sufficiently small. Since only small 
deformations are considered in the following, there is no distinction 
between the reference configuration and the current configuration. The 
displacement gradient H = Grad (u) ~ grad (u) is assumed to be ad- 
ditively decomposable into an elastic part H and a plastic part Hy. 
As a result, the additive decomposition also applies to the infinitesimal 
strain tensor 


e= sym(H) = sym(H.) + sym(H,) = Ee + ey. (3.6) 


Constitutive equations are required which describe the material-specific 
response of the body, e.g., to an external force. For material modeling 
it is crucial to describe the relation between stresses and strains so that 
the deformation of this material can be calculated. The general linear 
elastic anisotropic Hooke's law 


c = Cle.,] (3.7) 
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connects the stresses and strains by a linear relationship. The stiffness 
matrix C = Cine; & e; & ex ® ei is a fourth-order tensor with 81 com- 
ponents. Under consideration of hyperelasticity, the following symme- 
tries hold true. Due to the inherent symmetries of ø and e the stiffness 
matrix has two sub-symmetries C;;;; = Cjij; and Cijki = Cijik. Regard- 
ing the anisotropic elastic case, the symmetry Cijki = Cij;; and one of 
the sub-symmetries imply that only 21 elastic coefficients of the stiffness 
matrix are independent. This number of independent coefficients can 
be further reduced by material symmetries. 

Crystal plasticity frameworks are based on the assumption that plastic 
flow takes place by evolution of plastic slip on different slip systems. 
With the kinematical constitutive assumption 


H,= Y aM; (3.8) 
the evolution equation of the plastic strain e, results in 
N 
êp = $ ja M3, (3.9) 
a=1 


where the symmetric part of the Schmid tensor M5 := sym (Ma) is 
introduced. Based on the slip rate, the accumulated plastic slip on a 
slip system a is introduced as internal hardening variables by 


t 
Yaca | al dt. (3.10) 
0 


Finite deformation setting 


Based on the assumption of the existence of the plastic deformation 
p — P as the material isomorphism between two elastic materials 
(cf. Bertram (2008)) the deformation gradient F is multiplicative decom- 
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posed into the elastic deformation gradient Fe and the plastic deforma- 
tion gradient F', 
F = Grad (x) = F.F,. (3.11) 


The elastic deformation gradient F, = R.U. consists of the lattice 
stretch U, € Psym and the mean lattice rotation Re € Orth* of the 
crystal lattice. The plastic deformation gradient F, € Unim* repre- 
sents the distortion of the crystal lattice due to the formation and 
movement of dislocations. With the deformation gradient a material 
line element da = F dX, an area element da = det(F)F T dA and a 
volume element dv = det(F) dV can be transformed from the reference 
placement (upper case letters) to the current placement (lower case 
letters), respectively. 


Under the assumption of small elastic strains the St. Venant's law 


is introduced for the second Piola-Kirchhoff stress S.. Here, C denotes 
the elastic stiffness tensor in the intermediate configuration and 


E. = : (FIF. = I) (8.13) 


denotes the Green strain tensor. The relation of the second Piola- 
Kirchhoff stress in Eq. (3.12) and the Cauchy stress in Eq. (3.7) is given 
by o = det(F;')F.S.F. The multiplicative decomposition of the 
deformation gradient results in the velocity gradient 


L = grad (v) = FF} + FF FFF, (3.14) 


where the distortion rate tensor 


İp = FFF =Y uM (3.15) 


23 


3 Classical continuum mechanics theory 


consists of slip rates ^, on individual slip systems. In Eq. (3.15) the 
Schmid tensor for the slip system a 


Ma = da Q ña (3.16) 


was introduced in the lattice space. From this, a push forward to the 
Eulerian form can be performed 


M,-F.d,&F;ü,. (8.17) 


3.3 Mechanical and thermodynamical balance 
relations 


3.3.1 Motivation and mathematical preliminaries 


In the following it is assumed, as illustrated in Fig. 3.2, that a material 
volume Y is divided by a material singular surface S into two sub- 
volumes V+ and Y-, ie, V = Vt UV”. The outer normal vector of 
the surfaces 0V* and OV~ is denoted by n, while OV = 9V* UOY- 
holds true. The normal vector of the singular surface is denoted by 
ns and points from V- towards V^. The line of intersection of S and 
OV is denoted by ðS with outer normal vector m. The jump of an 
arbitrary quantity a at a singular surface is defined as Jar = at — a7, 
where at and a^ denotes the right-hand limit and the left-hand limit of 
a, respectively. 


Divergence theorem 


The divergence theorem can be used to transform a surface integral into 
a volume integral. For a material volume divided by a material singular 
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yr 


Figure 3.2: Material volume V, separated by a material singular surface S. 


suface it reads 
/ a-nda= | div(a) dos f Ja} ns da. (3.18) 
av y S 


Transport theorems 


The transport theorem specifies the time derivative of a volume integral 
of a bulk density Vy in a time-dependent volume Y = V(t). In the 
presence of a material singular surface it reads 


gf wave | EU w+ f Wyo nda | fev da, (3.19) 
dt Jy y ot av E 


where vg = vs: ns denotes the normal component of the velocity of 


the singular surface. 


The transport theorem for a surface density Vs in a time dependent 
material singular surface S = S(t) is given in Slattery (1990) as 


d 


Fl Vsdv = / (Ws + Vsdivs (v)) da, (3.20) 


where divs (-) denotes the surface divergence on S. For a detailed 
discussion on the surface divergence operator it is referred to Cermelli 
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et al. (2005), where the equivalence of 
divs (v) = divs (v!) — ksvt (3.21) 


is shown. Here, v! denotes the tangential component of the velocity v 
and ks = —divs (n) denotes the total mean curvature of the singular 
surface. Additionally, the surface gradient operator 


grad, (v) = Psgrad (v) , Ps := (I — ns 8 ns) (3.22) 
according to Cermelli et al. (2005) is introduced for a material surface. 
Equation of balance for volume and surface distributions 


The additive quantity V, which accounts for the bulk density Yy and 
the surface density Ws, is introduced 


v= | www | sav. (3.23) 
v S 


The rate of change of may occur due to three different causes: Flux, 
production and supply. The flux per element of area penetrated through 
OY and the flux per element of line penetrated through 0S are denoted 
by 9, and ®s, respectively. Additionally, the production densities 
pv and ps as well as the supply densities sy and ss, regarding the 
volume and the singular surface, are introduced. The general form 
of the equation of balance of a volume divided by a material singular 
surface reads 


S f vere [st | $y-ndo« | ®s-mdl 
dt Jy dt Js av as 


(3.24) 
+ | +) dw f (ps ss] da, 
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cf. Müller (1985). With the divergence theorem according to Eq. (3.18) 
and transport theorem according to Eq. (3.19) the local equation of bal- 
ance in regular points is obtained 


Vy + Vydiv (v) = py + sy + div ($,) s (3.25) 


Additionally, with the transport theorem for material singular surfaces 
according to Eq. (3.20) the local equation of balance in singular points 
results in 


Üs + V sdivs (v) — {Vy (v = vs)} "Ns 


. (3.26) 
= ps + ss + divs ($5) + {By} “Ns. 


For a material singular surface it follows that the velocity field of the 
singular surface is equal to the spatial velocity field of the body, i.e., 
v = vs and {uv} = 0. 


3.3.2 Balance equations 


Balance of mass 


It is assumed that the mass densities are conservative quantities, i.e., 
mass productions do not occur. Furthermore, it is assumed that mass 
fluxes across the exterior surfaces or lines are zero, and that mass can 
not be created inside the material by external actions, i.e., mass supplies 
are zero. The overview of all quantities of the balance of mass is given 
in Table 3.1, where o denotes the mass density of the bulk material and 
0s denotes the mass density of the singular surface. Consequently, the 
mass balance of a material body is given by 


d d 
= doas — — 0. . 
she v+ 5 f oda 0 (3.27) 
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Table 3.1: Overview of all quantities contained in the balance of mass. Balanced quantity 
V, production p, supply s and flux ® for bulk material and singular surface, respectively. 


With Eqs. (3.25) and (3.26) the local equations of balance of mass in 
regular points and singular points result to 


ò + odiv (v) — 0, (3.28) 


és + osdivs (v) = (3.29) 
Müller (1985) explains Eq. (3.29) by an illustration of curved imperme- 
able surfaces whose particles move perpendicular to the surface. In this 
study, however, the surface mass density os is assumed to be a constant 
on the dividing surface. Consequently, Eq. (3.29) simplifies to 


divs (v) — 0, (3.30) 
i.e., there is no local dilation of the surfaces, cf. Slattery (1990). 
Balance of linear momentum 


The balance of linear momentum is equivalent to Newton's second law 
of motion, which states that the rate of change of momentum of a body 
is equal to the forces acting on it. The applied forces are separated into 
body forces and contact forces. The body forces are denoted by b and 
bs for the body force per unit mass on the bulk material and surfaces, 
respectively. The contact forces are given with the Cauchy stress tensor 
c and the surface stress tensor os by the bulk material stress vector 
t = on and the surface stress vector ts = osm. The overview of all 
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Table 3.2: Overview of all quantities contained in the balance of linear momentum. 
Balanced quantity Y, production p, supply s and flux ® for bulk material and singular 
surface, respectively. 


quantities of the balance of linear momentum is given in Table 3.2. Con- 
sequently, the balance of linear momentum is given by 


d d 
gf ewt f omm 


= [ obau+ | osbsda+ f taa | tsdl. 
y S Oy Os 


By use of the balance of mass the local equations of balance of linear 


(3.31) 


momentum in regular points and singular points result to 


où = ob+div(o), (3.32) 
osù = osbs + divs (as) + 1o] ns. (3.33) 


Balance of angular momentum 


The balance of angular momentum states that when no external torques 
act on the material volume, no change of angular momentum will occur, 
i.e., the sum of all torques with respect to position £o acting on a body 
must be zero. The overview of all quantities of the balance of angular 
momentum is given Table 3.3. 
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Table 3.3: Overview of all quantities contained in the balance of angular momentum. 
Balanced quantity V, production p, supply s and flux ® for bulk material and singular 
surface, respectively. 


Et UE e d 
[zer erae fo eraa ferae 


zes [e axes [0] e a0) ests 


The balance of angular momentum reads 


d d 
a (ea) x ovdv+ 5 | (@- æo) x osvda 
-[ (æ — zo) x obav | (x — xo) x osbs da (3.34) 
v E 
+f (æ — a) x tda f (x — xo) x ts dl. 
av as 


With the assumption that all work on the body is the result of forces, 
not torques, it is shown in Slattery (1990) that the balance of angular 


momentum results in 
c —o!, geo. (8.35) 
Balance of kinetic energy 


The balance of kinetic energy is a corollary of the balance of linear 
momentum. The scalar multiplication of Eqs. (3.32) and (3.33) with the 
local velocity v results in the local equations 


Solo -v) = ob- v — o : grad (v) + div (c Tv) À (3.36) 
1 
zes: v) = esbs : v — os : grads (v) + divs (oiv )+Jo To]: ns. 


(3.37) 
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Table 3.4: Overview of all quantities contained in the balance of energy. Balanced quantity 
W, production p, supply s and flux ® for bulk material and singular surface, respectively. 


ee E 
rus [oure |ar| 
res] 


es anos | 0 | ases + osbs v 


From a comparison with the general forms of balance equations pre- 
sented in Eqs. (3.25) and (3.26) it follows that the kinetic energy is not a 
conservative quantity, since the production densities py = —o - grad (v) 
and ps = —os - grad, (v) are present. With the symmetry of stresses ac- 
cording to Eq. (3.35) and the definition of the infinitesimal strain tensor 
according to Eq. (3.6) it follows py = —e - é and ps =-os- (Psé). 


Balance of energy 


The energy of a system is given as the sum of internal energy and 
kinetic energy. The internal energy is denoted by e and es for the 
bulk material and singular surfaces, respectively. Being a conservative 
quantity, there is no production of energy. The supply of energy is 
given by heat supply densities ow and osws and supply densities of 
kinetic energy ob: v and osbs : v. The flux of energy is composed of 
non-convective heat transfer by the heat flux vectors —q and —qs and 
the contact energy transmissions ø ! v and ø 1v of the bulk material and 
the singular surfaces, respectively. The overview of all quantities of the 
balance of energy is given in Table 3.4. The balance of energy reads 


F eges) e x] mom 
T z e 3" v d es Les t 5 a 
= f ow+o-v) du+ | as (us +bs cv) da (3.38) 
y 


E 
+f (v-t—qen) do | (vts — qs: m) dl. 
ay Os 
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By use of the balance of mass the local equations of balance of energy 
in regular points and singular points result to 


" (: d jo v) = o (w +b- v) + div (ø 'w) - div (q) , (3.39) 


1 
os | ès + =ò- v | = os (ws + bs - v) + divs (odv) — divs (q 
s (és 5 ) s(ws+bs:v) s (osv) s (qs) (3.40) 


+ Jo w} -ns — 1q) ns. 
Balance of internal energy 


The local equations of balance of internal energy result by a subtrac- 
tion of the local equations of balance of kinetic energy from the local 
equations of balance of energy. Subtracting Eqs. (3.36) and (3.37) from 
Eqs. (3.39) and (3.40) results in 


ge = ow+a-é-div(q), (3.41) 
oss = OsWs t 0s: (Psé) — divs (qs) — Ja} “Ns. (3.42) 


Balance of entropy 


The entropy is not a conservative quantity, therefore, the entropy pro- 
duction densities op" and osp} are introduced for the bulk material 
and singular surfaces, respectively. For the entropy supply and the 
entropy flux the classical approach is chosen, cf., e.g., Coleman and Noll 
(1963), where the contribution to the balance of entropy results from the 
heat contribution divided by the absolute temperature 0. The absolute 
temperature is assumed to be continuous across singular surfaces, i.e., 
19} = 0. Consequently, the entropy density supplies and entropy flux 
vectors are given as listed in Table 3.5, where n and ns denotes the mass 
specific entropy of the bulk material and singular surfaces, respectively. 
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Table 3.5: Overview of all quantities contained in the balance of entropy. Balanced 
quantity V, production p, supply s and flux ® for bulk material and singular surface, 
respectively. 


= [sor aw+ f as (2 +08) da (3.43) 


By use of the balance of mass the local equations of balance of entropy 
in regular points and singular points result to 


2000 Bu suc PT 
on = es + op" div (4) , (3.44) 
1 : 1 
0sl)s = 05 + 0sps — divs (38) -= ziel "Fs. (3.45) 


3.3.3 Clausius-Duhem inequality 


The second law of thermodynamics states that, in a closed system, the 
minimum time rate change of entropy is equal to the rate of entropy 
transmission. Consequently, the total production of entropy is non- 
negative. In Coleman and Noll (1963) it is stated that it is necessary and 
sufficient that the specific entropy production is non-negative at all ma- 
terial points and at all times, i.e., op" > 0 and osp% > 0. With D :— op"0 
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the dissipation of bulk material is defined, whereas Ds :— osp20 de- 
fines dissipation of singular surfaces. From a Legendre transformation 
it follows 


V —e-—m6, (3.46) 
Us — es — ns? , (3.47) 


cf., e.g., Levine (2009), where the specific Helmholtz free energy of the 
bulk material y and the specific Helmholtz free energy of singular sur- 
faces ws are introduced. With Eqs. (3.41), (3.42), (3.44) and (3.45) the 
resulting inequalities, called Clausius-Duhem inequalities for regular 
points and singular points, read 


D --—wo)-— mn +0 -è Ža- grad (0) > 0, (3.48) 


, : . 1 
Ds = Qs s osnsé t Os: (P sé) — 945 $ grads (0) > 0. (3.49) 


The mass density of the bulk material, the mass density of the singular 
surfaces, as well as the temperature, are per definition non-negative, i.e., 
o > 0,0s > 0 and 0 > 0. 


3.4 Summary and conclusion 


The deformation of a body is influenced by both external loads and and 
internal loads. Balance relations are introduced which provide state- 
ments about the connection between the rate of change of a mechan- 
ical or thermomechanical quantity based on production, supply and 
flux contributions. In continuum mechanics the mechanical balance 
relations of mass, linear momentum, angular momentum and kinetic 
energy are considered. With the balance of energy, internal energy and 
entropy the Clausius-Duhem inequalities are obtained by means of the 
second law of thermodynamics. Hereby, regular points and singular 
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surfaces are considered in order to model the mechanical behavior of 
the bulk material as well as the mechanical behavior of GBs in polycrys- 
talline materials. This provides the basis to derive thermodynamically 
consistent flow rules for the plastic deformation inside the bulk mate- 
rial and at the GBs. In order to account for the nonlocal mechanical 
behavior of polycrystalline materials, which results from interactions 
between dislocations and GBs, however, additional DOF have to be in- 
troduced in the framework. In the following chapter a gradient crystal 
plasticity theory is established in the context of small deformations by 
considering the plastic slips as additional DOF. 
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Chapter 4 


Gradient crystal plasticity theory 


4.1 Introduction and motivation 


For polycrystalline structures classical continuum theories fail to repro- 
duce the experimentally observed influence of grain size on the yield 
stress, cf., e.g., Armstrong et al. (1962); Conrad et al. (1967). In order to 
account for work-hardening in polycrystals, Estrin et al. (1998) divides 
the cell structure into cell walls and cell interiors and formulates the 
resolved shear stress as a function of two dislocation densities, respec- 
tively. At the core of this and most phenomenological plasticity models 
are the evolution equations for the dislocation densities. A commonly 
used evolution law for FCC crystals accounts for dislocation storage 
and dislocation annihilation via the Kocks-Mecking relation, cf. Meck- 
ing and Lücke (1970); Kocks and Mecking (2003). While theories de- 
pendent on randomly distributed SSD are able to predict the overall 
mechanical behavior of homogeneous materials, plastically inhomoge- 
neous deformations on the microscale require special consideration. To 
consider the work-hardening induced by interactions between disloca- 
tions and GBs, Ashby (1970) introduces the density of GND which is 
assumed to scale proportionally with the amount of plastic slip on each 
slip system. Based on this approach various strain gradient plasticity 
theories have been developed. In the early work of Aifantis (1987) the 
macroscopic back stress is formulated with regard to the Laplacian of 
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the plastic slip. More recent works include experimental investigations 
(e.g., by Fleck et al. (1994)), applications to polycrystals (e.g., by Ohno 
et al. (2008)), finite deformations (e.g., by Gurtin (2008)) and dissipative 
solids (e.g., by Miehe (2011)). 

By introducing generalized stresses power conjugate to the in-plane 
gradients of plastic slip, Gurtin (2002) develops a thermodynamically 
consistent framework based on the principle of virtual power. Follow- 
ing the same approach, a thermomechanically coupled gradient theory 
is derived in Anand et al. (2015). While Gurtin (2002) uses an exten- 
sion of the principle of virtual power to derive the governing equations 
for the gradient-enhanced plasticity model, Svendsen (2002) formulates 
the constitutive models by using an invariance consideration of the 
extended balance of total energy. In Prahs and Bóhlke (2020b) the in- 
variance of an extended energy balance is considered with respect to a 
change of observer. A contribution of kinetic energy associated with the 
scalar-valued slip is taken into account. This chapter gives an overview 
of the required assumptions in order to get a consistency of the model 
approaches by Gurtin (2002) and Prahs and Bóhlke (2020b). 


4.2 Model approach based on extended bal- 
ance equations 


4.2.1 Extended balance equations 


Extended energy balance 


For each slip system o the plastic slip ya is introduced as an additional 
DOF. Each additional scalar-valued DOF extends the balance of energy, 
given in Eq. (3.38), by an additional kinetic energy density, an addi- 
tional volume specific mechanical power and an additional surface spe- 
cific mechanical power. The spatial velocity of the additional DOF is 


38 


4.2 Model approach based on extended balance equations 


Table 4.1: Overview of all quantities contained in the extended balance of energy. 
Balanced quantity v, production p, supply s and flux ® for bulk material and GBs, 
respectively. 


vm In 0Ayaya/2 a j D 
oses + os(v- v)/2 0sus--0sbs:v | -qs + 08V 
at: ou 05 Ásza 34/2 TOS pat bsyaYSa VY > EsaYsa 


given as the time derivative of plastic slip, i.e., Ya. For each slip system 


a the micro-inertia force A,., the generalized micro-body force bya and 
the generalized micro-stress €, are introduced, cf. Prahs and Böhlke 
(2020b). Additional contributions regarding the GBs are denoted by 
Asya, bsya and £5,, respectively. An overview of all contributions to 
the balance of energy extended by plastic slips is given in Table 4.1. 
With Eq. (3.24) the extended balance of energy is given by 


al + dy) d +f yee. a 
di nd € p^ v U dt OS | €s 2 a 
d fix " ad Pi T 
Hg] otn gy [5 asado 


= [oet au+ | as (ws +bs v) da 
y S 


N N 
«f os; booa au+ | es Y bs asa da 
Y a=l S a=1 
+f reed (v: ts — qs: m) dl 
OV in 


N 
+f Se n) Vex da «f. > (Esq: m )^se dl. 
OV a= 


(4.1) 
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In Prahs and Bóhlke (2020b) the local form of the balance of mass and 
the local form of the balance of linear momentum are obtained for reg- 
ular points under the assumption of invariance of the extended energy 
balance with respect to a change of observer. The extension to a material 
including singular surfaces is given in Prahs and Bóhlke (2020a). From 
the invariance considerations it follows that neither the local balance 
of mass nor the local balance of linear momentum is affected by the 
extension of the balance of energy by additional scalar-valued DOF. By 
use of the balance of mass, given in Eq. (3.28), the local equation of the 
extended balance of energy in regular points results to 


1 iz bes 
[epe eX e Yet) 


a=1 a=l 


N 
=o (+ +b:v+), s (4.2) 


a=l 


+div (a! v) — div (q) + > div (Eaa) 


Analogously, with Eq. (3.29), the local equation of the extended balance 
of energy in singular points follows 


1 ig D 
OS (5 + 3? -UF 3 5 ASya Sa} Sa + 3 X Ast] 


a=1 a=1 
N 
= 0s c +bs:v+ 5 iss 
a=1 
" (4.3) 
+divs (og v) — divs (qs) + >> divs (&saYsa) 
a=1 
N 
He 'vļ} -ns — dab ns ++ 5 [eater ns. 
a=1 
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Extended internal energy balance 


By subtracting the local form of the kinetic energy, given in Eq. (3.36), 
from the extended balance of energy, given in Eq. (4.2), the local form 
of the extended internal energy in regular points is obtained 


N 
1 : 
oÈ — Q0 (+ us 5 c m 2 (Aaa E A) ) io 

BS " (4.4) 


to 6 - V div (€,4a) - div (q) . 


a=1 


Analogously, the local form of the extended internal energy in singular 
points is obtained by subtracting Eq. (3.37) from Eq. (4.3) 


N 
1 
ests = 0s (s + > (iss. 75 (Asya¥sa + Asaisa)) 3 
a=l 
N 
tos: (Pse) +) divs (&sa/sa) - divs (qs) (45) 
a=l 


+ > 1£. 3] ns -Iq} ns. 


Simplifying assumptions 


Following the work of Prahs and Bóhlke (2020a;b), effects associated 
with micro-body forces bya and bsya and micro-inertia Aya and Asya 
are neglected. Consequently, it is assumed that there are no body forces 
or kinetic contributions associated with the plastic slips. Additionally, 
the surface stress vector of the GBs ts, i.e., the contact force per unit 
length at the curve S, as well as the additional contribution associated 
with the plastic slip £5, : m, are neglected as well. This simplification 
is motivated by the neglection of GB sliding within this work. 
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4.2.2 Clausius-Duhem inequality 


Under consideration of these assumptions the Clausius-Duhem in- 
equality for regular points and singular points result to 


N 
: : 1 
D = —op — od +0 -è+ ) div (Eaa) — gd grad (0) > 0, (4.6) 


a=1 


N 
Ds = -0stbs — 05150 + 5 {E40} ns - ads grad, (0) >0. (4.7) 


a=l 


In order to derive the reduced dissipation inequalities constitutive as- 
sumptions are required for the free energy of the bulk material and the 
free energy of the GBs. 


Bulk material 
The Helmholtz free energy of the bulk material is assumed in the form 
Y = SP Ee Yaca» grad (Ya) ,8) . (4.8) 


With this, the derivation of the Helmholtz free energy is given by 


_ ab. æ 
v Oe "Ee Oe P 
N (4.9) 
op Ow Ow . 
+ 5 = Yaca + EIA A grad 6) + 590 


Hereby, slip fields are varied independently since they represent addi- 
tional fields of DOF. The derivation with respect to the accumulated 
plastic strain can be further specified by 


av. dv 
Maca ea B 


al = x sen (^) Ya - (4.10) 
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Inserting Eq. (4.9) in the Clausius-Duhem inequality of the bulk mate- 
rial, given in Eq. (4.6), results in 


_ do) Š dow 
dal ie RAC >) mtt 


: Op : Op s |a 
+) (div ie p. EP (Ya) — de, L3 7& (4.11) 
a=1 
dow\, 1 
T (on + 2e) Ó — 2a: grad (0) > 0, 


Ve, é, 0, Å, Yas Yas grad (Ya) ‚grad (Ya) $ 


where the evolution equations of the plastic strains, Eq. (3.9), is taken 
into account. Assuming that the Cauchy stress ø as well as the gradient 
stresses €, are purely energetic, the potential relations are obtained 


| dow 
o = 9m (4.12) 
_ 0p 
Sa = Gerad (34) ' Va € {1,...,N}, (4.13) 
0 
on = me . (4.14) 


This leads to the reduced dissipation equation of the bulk material 


—— 
thermal contribution 

(4.15) 
which can be decomposed in a mechanical contribution and a thermal 


mechanical contribution 


contribution. Here, both contributions of the dissipation are considered 
to be uncoupled. With Fourier's law for isotropy 


q — —Aegrad (0) , (4.16) 
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cf., e.g., Bergman et al. (2017), where Ae denotes the thermal conduc- 
tivity, the positivity of the thermal contribution to the dissipation is 
ensured. Consequently, the reduced dissipation equation of the bulk 


material is given as 


N 
. o : 0 B es 
Dred = p» (div (£4) — x sen (Ya) — v Mm) Ja 20. (4.17) 


a=l 


Grain boundaries 


The Helmholtz free energy of the GBs is assumed to depend on the 
temperature as well as on the right-hand limit and the left-hand limit of 
the plastic slips, i.e., 

bs = Vs(Ya Ta 0). (4.18) 


Slip fields are varied independently which results in 
00 vf) OU Os ; 
ds =~ (Sei + ei) +6. (4.19) 


Oye ^ a 


a=l 


With Eq. (4.7) this results in the dissipation of the GBs 


N p 
Ds =). (et ns - 25888): zx ns - Pests) 4, 


a=l 


- (asus 25578) — fas grads (0) > 0, 


V0,0, y AT yn g 
(4.20) 


From this the potential relation 


esns = ———— (4.21) 
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is obtained. In analogy to Eq. (4.16) the Fourier's law for GBs 
ds = -Aosgrad; (0) , (4.22) 


with the thermal conductivity of GBs Aes, ensures the positivity of 
the thermal contribution of the dissipation. Consequently, the reduced 
dissipation equation of the GBs is given as 


Y dos 
Dsrea = X ( aig js ar 
g= m (4.23) 


N 
dos = 
+5 (-& ns - Ss) 5, 2 0. 


4.3 Model approach based on principle of vir- 
tual power 


4.3.1 Principle of virtual power 


In the following, based on the principle of virtual power, so called mi- 
croforce balance equations are obtained for each slip system. It is Shown 
that the microforce balance equations are equivalent to nonlocal yield 
conditions for the individual slip systems. 

For the evolution of the polycrystalline body, composed of bulk mate- 
rial and GBs, a microscopic and a macroscopic force system is consid- 
ered. The microscopic system in Cermelli and Gurtin (2002) is defined 
by contributions of both, bulk material and GBs. The bulk material 
contributions are a lattice stress ø that expends power over the lattice 
strain-rate € and, for each slip system o, a scalar-valued microstress 
Ta that expends power over the slip-rate Ja and a vector-valued mi- 
crostress £, that expends power over the slip-rate gradient grad (ya). 
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Regarding the GBs, microtractions II} and II, expend power over the 
right-hand limit and left-hand limit of the slip rate yt and yz, respec- 
tively. In Cermelli and Gurtin (2002) an additional GB contribution is 
given by a macroscopic internal force 7 conjugate to the rate of dis- 
placement discontinuity {ù}. Cavitation at the GB is excluded, i.e., 
fu} ns = 0, however, the jump {u}, which is tangential to S, intro- 
duces GB slip. This effect is neglected in this work. Based on the 
individual contributions the internal power is given as 


N 
Pis x c6 Y (Tata + E, grad (Ja) | de 
v a=1 


(4.24) 


N 
«f XO tat 154) da. 
S a=1 


The power of external forces is constructed based on the work of Gurtin 
(2002). The macroscopic system is defined by a body force density 
ob and a macroscopic surface traction t, both associated with à, and 
further, for each slip system o, by a microscopic surface traction E, 
associated with the plastic slip rate 4. Consequently, the power of 
external forces is given as 


N 
Pu = | ob-wav | Pado f 2 ajada, (4.25) 
Y 8Y, 8Ys 4-1 


where OY; denotes the Neumann boundary corresponding to t and ÖV= 
denotes the Neumann boundary corresponding to =. The principle of 
virtual power states that given any generalized virtual velocity óv and 
any subregion the corresponding internal and external virtual powers 
are balanced. For Eqs. (4.24) and (4.25) the generalized virtual velocity 
is introduced as the list dv = [öü, de, 0*y, }, where dt, de and ya denote 
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the virtual counterparts of à, & and ^4. The power balance reads 


N 
f G det Y (adja + Eq: grad un du 


a=1 


N 
+ J XO (m4 +My ey) da (4.26) 
S a=1 


ob: sid | t- sada- | E05, da = 0. 
=, OY. e» 


With the balance of angular momentum, i.e., o = 0 T the application of 
the divergence theorem, given in Eq. (3.18), results in 


f Lo a — div (£,)) Ya w- f (ob + div (o)) - ù dv 
+ h Y (TA, + HZ —1£,oy4]- ns) da — (427) 


«f > - Eq ) ida | (en - t) Suda — 0. 
OYz OVt 


Eq. (4.27) has to be satisfied for all (04, ôy} and all subbodies, i.e., for 
bulk material as well as for GBs, separately. This results in the following 
local field equations, boundary conditions and GB conditions 


Ta — div(£,) =0, Va € {1,...,N}, VeeV, (4.28) 

gb + div (o) =0, Vrcy, (4.29) 
Zt- Et- ng =0, Va € {1,...,N}, VaeS, (4.30) 
=, téz ns =0, Va € {1,..., N}, Vaes, (4.31) 
£,:n— 8, — 0, Va € {1,...,N}, VeeöVs, (4.32) 
on-t=0, Va c 0V. (4.33) 
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With the pill-box theorem of Müller (1985) ns = nj; = —n,, holds true 
for the singular surface normal vectors of the right-hand subregion V+ 
and the left-hand subregion V^. Consequently, Eqs. (4.30) and (4.31) 
can be written as 


= 
= 


+ _¢t nt =0, Va € {1,..., N}, WreS. (4.34) 


4.3.2 Clausius-Duhem inequality 


In the following, the Clausius-Duhem inequality is derived based on 
the work of Gurtin (2008), where a purely mechanical framework is con- 
sidered. Subsequently, a comparison to the thermomechanical frame- 
work based on extended balance equations is made. The derivation of 
the dissipation inequality is based on a second law in which, for any 
subregion of the body, the temporal increase in free energy is less than 
or equal to the power expended on it. Thus, for the bulk material and 
the GBs it follows 


/ Bi dus Pu: / osbsda< Perl), (435) 
y T 


where Pex:(V) denotes the power expended on the bulk material and 
Pext(S) denotes the power expended on the GBs. With the power bal- 
ance and Eq. (4.24) the dissipation of the bulk material D and the dissi- 
pation of the GBs Ds are obtained by localization 


N 
D= -oh +o è+ Y (Tata tq - grad (34)) = 0, (4.36) 
a=1 
] N 
Ds = —osis +), (Id +54) 2 0. (4.37) 


a=1 
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In order to derive the reduced dissipation inequalities constitutive as- 
sumptions are required for the free energy of the bulk material and the 
free energy of the GBs. 


Bulk material 


In the following, the Helmholtz free energy of the bulk material is as- 
sumed in the form 


V = V(Ee, Yaca, grad (Ya)) - (4.38) 


Inserting in the Clausius-Duhem inequality results in the bulk material 
dissipation inequality 


N 
den? | (439) 
D» En - = Ga) Ma) ia 20, 


V ££, Yas Yas grad (Ya) ‚grad Ya) « 


Assuming that the Cauchy stress ø as well as the gradient stresses €, 
are purely energetic, the following potential relations are obtained 


| dow 
ga (4.40) 
107020) 


Le 1,...,N}. 441 
gad Ga) Va € {1,..., N} (4.41) 


E47 


This leads to the reduced dissipation equation of the bulk material 


N 
a 0. 8 
Duae (m = sot sgn (44) - a TE 420. (442) 
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Grain boundaries 


The Helmholtz free energy of the GBs is assumed to depend on the 
right-hand limit and the left-hand limit of the plastic slips, i.e., 


Vs = slyt, va). (4.43) 


Slip fields are varied independently which results in the dissipation of 
the GBs 


D ler Doss) 4 (= 00sUs EBEN" 
Sed = >, zz dt Ea = -a Ya Z0. (4.44) 


a=1 a=l 


4.4 Comparison and simplifications 


To compare the model approaches the following assumptions are made 


* Helmholtz free energy temperature independent 


v z v(0), Vs # Vs(0) (4.45) 


* Isothermal processes 
6=0 (4.46) 


* Homogeneous distribution of temperature 


grad (0) — 0, div (q) 2 0 (4.47) 


* No heat sources 
w=0, ws =0 (4.48) 


With Eq. (4.45) and the potential relations of Eqs. (4.14) and (4.21), the 
Legendre transformations of Eqs. (3.46) and (3.47) result in Y% = e and 
Ws = es for the bulk material and GBs, respectively. Consequently, 
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with the made assumptions, the balance of internal energy, given in 
Eqs. (4.4) and (4.5), results in 


N 
oè =0- è+ V div (faja) , (4.49) 
" a=1 
Osés = 5 1&5]: ns. (4.50) 
a=1 


With Eqs. (4.6) and (4.7) it directly follows that both, the bulk material 
and GBs, are dissipation free, i.e., 


N 
D --g)ceo-6- V div (Eata) - 0, (4.51) 
a=1 
. N 
Ds = -eshs +X Era} ns — 0. (4.52) 
a=1 


An in-depth discussion of the vanishing dissipation of a purely me- 
chanical theory within the framework of continuum thermodynamics 
is given in Prahs and Bohlke (2022) for extended continua with addi- 
tional non-thermal DOFs. Whereas in Prahs and Bohlke (2022) singular 
surfaces are not considered, it is shown by Eq. (4.52), that with the as- 
sumptions of Eqs. (4.45)-(4.48) the dissipation on the singular surfaces 
vanishes. The vanishing dissipation on singular points is also obtained 
in Prahs and Bohlke (2020a). 

The comparison of the reduced dissipation inequalities obtained by 
extended balance equations, given in Eqs. (4.17) and (4.23), with the 
reduced dissipation inequalities obtained by the principle of virtual 
power, given in Eqs. (4.42) and (4.44), results in m. = div (€,) and 


Fa =a ny- These are the local field equations obtained from the 


=i 


principle of virtual power, i.e., Eqs. (4.28) and (4.13), which are con- 


tained implicitly in the model approach of extended balance equations. 
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4.5 Bulk material flow rule 


In the following, a bulk material flow rule is constructed which ful- 
fills the reduced dissipation inequality of the bulk material, given in 
Eq. (4.6). With the additive decomposition of strain into an elastic 
and plastic part, cf. Eq. (3.6), and the general linear elastic anisotropic 
Hooke's law, cf. Eq. (3.7), the potential relation 


Op _ 


— = 4. 
Be, c (4.53) 


is obtained from Eq. (4.12). Furthermore, the following notation is intro- 
duced for the derivation of the free energy density by the accumulated 
plastic slip of slip system a 


o 
Ba = oy 


= 4.54 
Maca ( 3 ) 


With the result, that the bulk material is dissipation free, cf. Eq. (4.51), 
the reduced dissipation inequality of the bulk material, given in Eq. (4.6), 
results in 


(Ta + div (£a )) sgn (Ya) — Ba) Yaca =0, Va € (hbe Nye (4.55) 


Here, Eq. (3.10) is taken into account stating that 44.4 = |Y«|. Assuming 
sen (Ya) = sgn (Ta + div (€,)) leads to 


(Ira + div (E)| — Ba) Vaca =0, Vo € {1,..., N}. (4.56) 


To characterize plastic flow a yield criterion function Ya < 0 is defined 
for each slip system. For details on yield criterion functions it is referred 
to Simo and Hughes (1998). In this work, the rate-independent yield 
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criterion of slip system a is specified by 
Pa :— |Ta + div (€,)| — Bo. (4.57) 


Here, ba characterizes the isotropic hardening of the yield surface and 
div (£,) the kinematic hardening defining the center of the yield surface. 
Eq. (4.56) then reads 


Payaca — 0, Vaetl,...,N}. (4.58) 


From this the loading and unloading conditions can be formulated for 
each slip system a in the Kuhn-Tucker form with 


PaYaca = 0, Pa < 0, Vaca = 0. (4.59) 


In contrast to viscoplasticity in the multi-slip case an active set search 
algorithm has to be used in the rate-independent case. The set of active 
slip systems is generally non-unique. For a comparison of algorithmic 
approaches it is referred to Miehe and Schröder (2001). 


4.6 Summary and conclusion 


In this work, the modeling of plastic slip under consideration of dis- 
location pile-up at GBs is realized by an extended energy balance ac- 
cording to Prahs and Bóhlke (2020a;b). The framework is compared to 
the model approach based on an extended principle of virtual power 
according to Gurtin (2002). The two model approaches are equivalent if 
isothermal processes, a homogeneous distribution of temperature and 
a temperature-independent free energy for the bulk material and the 
GBs are assumed. It is shown, that under these assumptions the bulk 
material and the GBs are dissipation free. The following chapter is 
dedicated to the investigation of the GB flow rule. 


53 


Chapter 5 


Comparison of grain boundary 
models 


5.1 Introduction and motivation 


In the early works of gradient plasticity theories, GBs act as insur- 
mountable obstacles for plastic slip, cf., e.g., Gurtin (2002). Gurtin 
and Needleman (2005) discuss microhard and microfree boundary 
conditions within a gradient crystal plasticity framework. A microhard 
GB state corresponds to the situations when GBs act as barriers to 
dislocation motion, hence the plastic strain within the GBs vanishes. 
The microfree GB condition implies that the GBs allow free dislocation 
transmission so that the microtractions are prescribed to be zero, i.e., 
€,: ns =0 Va. In order to model a GB behavior in between those 
limits, slip transmission criteria have to be introduced. The plastically 
active state of the GBs is referred to as microcontrolled in the follow- 
ing in order to highlight the controlled evolution of the GB slip by 
additional field equations. 

A lot of experimental work has been done in order to detect GB resis- 
tance to slip transfer, e.g., by the use of nanoindentation. In Soer et al. 
(2005) a displacement jump is observed during the loading part of the 
indentation. It is shown in Wang and Ngan (2004) that this so called 
pop-in effect has a close relationship with the misorientation across the 
GB. The observed trend that the slip transmission ability of a GB is 
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mainly governed by the misorientation between slip planes and slip 
directions across the GB is supported by Wo and Ngan (2004). The 
direct observation of the interaction of individual dislocations with GBs 
of Kondo et al. (2016), however, shows that beside the geometric effects 
of misaligned slip systems, local structural stabilization effects cause 
the dislocation impediment at GBs. 

It is stated in Wo and Ngan (2004) that nanoindentation are not an 
effective method to probe the slip transfer properties of GBs. Especially 
in FCC materials, where there is little dislocation locking and little re- 
sistance to slip transfer across GBs, a GB pop-in is not observed, cf. 
Britton et al. (2009). The efficient mechanical characterization of GBs 
using nanoindentation is therefore mainly restricted to BCC metals. For 
the characterization of GBs within BCC metals the mechanical response 
to nanoindentation is recorded as a function of the distance to the GB, 
cf., e.g., Soer and De Hosson (2005); Ohmura and Tsuzaki (2007). The 
characterization is strongly limited by the instrumented indentation 
and the recorded resolution. More recent work, e.g., Kalidindi and 
Vachhani (2014), therefore address these problems and obtain advances 
by automation of measurement and data analysis. Besides nanoindenta- 
tion, in-situ electron microscope measurements and diffraction-contrast 
electron tomography can be used to investigate dislocation work hard- 
ening effects, cf., e.g., Lagow et al. (2001), or dislocation-GB interac- 
tions, cf., Kacher and Robertson (2014). Interactions of dislocations at 
GBs, however, are complex with multiple dislocation systems emitted 
simultaneously. Therefore, models and simulations of slip transmission 
across GBs are essential in order to understand the slip system activa- 
tion within polycrystalline materials. 

Various realizations of slip transmission criteria in simulations of metal 
materials exist. In Hamid et al. (2017) a GB model within a continuum 
dislocation dynamics framework is constructed. For each slip system 
pair a slip transmissivity number is computed based on the angle be- 
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tween the slip planes and the angle between the slip directions. This 
geometrical condition was introduced in Werner and Prantl (1990). The 
GB strength is assumed as a linear function of the slip transmissivity 
number. Alipour et al. (2020) construct a GB model based on the same 
geometrical condition of Werner and Prantl (1990) within a gradient- 
extended framework. Hereby, an overall measure of all slip system 
combinations is constructed. In contrast to Hamid et al. (2017) the as- 
sumed form of the GB yield strength is able to predict the limiting case 
of a microhard GB. Other simulation approaches for slip transmission 
phenomena are dislocation-density based crystalline plasticity formu- 
lations, cf., e.g., Shi and Zikry (2011), the modification of a classical 
crystal plasticity model where slip transfer is permitted only for low- 
angle boundaries, cf., Haouala et al. (2020), or atomistic simulations, cf. 
Spearot and Sangid (2014). The geometrical slip transmission criteria 
have in common that transmission prevails at low-angle GBs. It is 
shown in Liu et al. (2020) that in contrast to this, at high-angle GBs 
dislocation source activation occurs. This effect, however, is not consid- 
ered in this work. For a detailed overview on slip transmission criteria 
in metals it is referred to the review article of Bayerschen et al. (2016). 

In the following, the GB interface energy of Gurtin (2008), dependent 
on Nye's dislocation density tensor, is investigated. The surface energy 
incorporates both the orientation of the GBs relative to the crystal lattice 
of the adjacent grains and the misorientation of slip systems between 
adjacent grains. Intra-grain interactions and inter-grain interactions of 
dislocations on different slip systems provide a physically meaning- 
ful framework for polycrystalline structures. The numerical simula- 
tion of Ozdemir and Yalcinkaya (2014) for the two dimensional case, 
or Gottschalk et al. (2016) and McBride et al. (2015) in three dimensions 
provide valuable insight into the GB model of Gurtin (2008). The plastic 
slip gradient at the GB is only dependent on the respective energy of 
the GB. The plastic slip rates at the GB are further assumed to be inde- 
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pendent. As shown in this work, for the limiting case of a vanishing 
grain misorientation, this only results in a single-crystal response for 
equally sized grains. In order to overcome this shortcoming, a second 
approach to derive the GB condition is introduced. Peng et al. (2019) 
couple the plastic slips for the incoming slip systems with the outgoing 
slip systems by kinematic relations. Hereby, however, a continuous 
distribution of plastic slip across GBs is assumed. In this work the 
general case of discontinuous plastic slip distributions is considered. 
An additional GB flow rule motivated by Wulfinghoff et al. (2013) and 
Erdle and Bóhlke (2017) is introduced based on the gradient of plastic 
slip. Instead of the coupling of slip gradients of adjacent grains with 
a constant material parameter the coupling is formulated based on a 
misorientation function. 

For the analytical investigation of the GB modeling approaches the 
framework is applied to the shearing of a three-phase periodic lami- 
nate structure for the special case of single slip. A similar problem, 
Le. shearing of a two-phase structure consisting of an elasto-plastic 
phase and an elastic phase, has been investigated by Sedlácek and 
Forest (2000) for a dislocation-based model of plasticity. Forest and 
Sedláček (2003) consider the identical boundary value problem in order 
to compare dislocation-based, Cosserat and strain-gradient models of 
crystal plasticity. In Forest (2008) the extension to a cyclic deformation 
is made in order to investigate kinematic hardening within the models. 
Additional applications to the two-phase laminate have been made 
in Cordero et al. (2010); Aslan et al. (2011) for the investigation of 
size-effects and in Forest and Guéninchault (2013); Wulfinghoff et al. 
(2015) for the investigation of free energy potentials in gradient crystal 
plasticity. Erdle and Bóhlke (2017) consider a two-phase laminate 
consisting of two elasto-plastic phases in order to investigate GB effects 
within a gradient crystal plasticity framework. Prahs and Bóhlke 
(2020a) reintroduce a purely elastic third phase for the investigation 
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of the influence of energy parameters with energetic GBs according 
to Gurtin (2008). Hereby, the special case of coherent slip systems, 
orthogonal to the GBs is considered. In this work the analytical solution 
of the three-phase laminate structure is constructed for a variable 
misorientation between the two adjacent elasto-plastic grains. The 
effect of the misorientation is investigated with special interest of a 
single crystal consistency check for coinciding slip systems. 


5.2 Grain boundary flow rule 


5.2.1 Grain boundary energy 


In order to derive the GB conditions that fulfill Eq. (4.23) an expression 
of the specific free energy of the GB vs is required. In this work a GB 
energy is chosen based on a measurement of GND. An experimental in- 
vestigation of GND distribution in an aluminum bicrystal is performed 
in Sun et al. (1998; 2000). Here, it is shown that at first GND concen- 
trate near the interface and redistribute with higher strain levels. This 
motivates a GB flow rule based on the amount of GND concentration 
at the GB and the activation of GB slip at higher strain levels. In order 
to characterize GND, Cermelli and Gurtin (2001) introduce the geomet- 
ric dislocation tensor G = Fcurl (Fp) /det(F',). G measures the local 
Burgers vector, which represents the closure deficit of circuits deformed 
from a perfect lattice, per unit area in the microstructural configuration 
within a finite deformation setting. A gradient theory with a bulk mate- 
rial defect energy dependent on G is presented in Gurtin (2002). Here, 
G = Y, 0da & da + o, (no x da) & dq is expressed in terms of pure 
screw dislocations o? and pure screw dislocations o^. Additionally, for 
a small deformation theory, G = curl (Hp) is given. A small deforma- 
tion theory based on this is presented in Cermelli and Gurtin (2002) 
where the specific defect energy v = c || G || ? /2is chosen ina quadratic 
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isotropic form with material constant c. Based on this defect energy 
approach Gurtin and Needleman (2005) construct a possible GB condi- 
tion. The Burgers vector density of the GB, measured per unit length, is 
characterized by the GB Burgers tensor Gs = (H | (nsx), which reads 
in index notation (Gs)i = {(Hp)ijfejx(ns)x. Consequently, a GB is 
referred to as defect-free if | HE (ns x) = 0 and the Burgers vector pro- 
duction within the GB is given as Gs = || {H,}(nsx) ||. Gurtin (2008) 
specifies the GB energy in a quadratic form 


1 
ests = 55| Gs IP , (5.1) 


based on the GB Burgers tensor. Here, the GB strength « is introduced 
as a material parameter. With Eq. (3.8) the magnitude of the GB Burgers 
tensor can be formulated in terms of slip-interaction moduli of each slip 
system pair a, B as 


N N 
2 t 
Gs? =X >) CH + Cog tare -20i vere, 62 
a=1 B=1 


cf. Gurtin (2008). The slip-interaction moduli 


ott = (ai di) (nz x ns) ; (nj x ns) (5.3) 
Cas = (a; dj) (n; x ns) : (n; x ns) ; (5.4) 
Cis = (at : d;) (nz x ns) : (n; x ns) (5.5) 


depend on the Schmid tensors Ma = da & n, of the two grains and the 
GB normal vector ns. Cy and Cag can be interpreted as intra-grain 
interaction moduli for the respective grain which represent interactions 
between slip systems within the respective grain. Cis is interpreted 
as an inter-grain interaction modulus which represents interactions be- 
tween the two grains. 
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5.2 Grain boundary flow rule 


In the following dissipation free GBs are considered, cf. Eq. (4.52). With 
a discontinuous slip distribution across GBs the reduced dissipation 
inequality of the GBs, given in Eq. (4.23), results in N equations for 
the 2N GB slip values. Two approaches are presented in order to solve 
this underdetermined system of equations. 


5.2.2 Approach I: Independent plastic slips 


One approach, here called approach I, is to assume independent slip 
rates at the GB. With this assumption Eq. (4.52) can directly be solved to 
the 2N unknown GB slip values. With the assumption of independent 
slip rates at the GB it follows Approach I (independent plastic slips) 


a y "WE 
gt ns=es 22S &(chwi-0nas). 66 
Ya 8-1 
- ð OE 
“ta ns me ra). 6D 


B=1 


5.2.3 Approach II: Coupled plastic slips 


In contrast to approach I, in the following the plastic slip rates at the GB 
are not considered to be independent. An additional flow condition 
based on a misorientation function f (Qi; ‚ns) is introduced. QU 
describes the misorientation between the slip systems a and 5 of the 
adjacent grains, e.g., by three Euler angles. In this work it is proposed 
that the misorientation function correlates to the amount of dislocations 
that piled-up at the GB, i.e., to the discontinuity of the plastic slip gra- 
dient at the GB. If two slip systems a and 3 hold no misorientation, 

n — I holds true and no dislocations pile-up at the GB, which results 


in a continuous gradient of slip, i.e. (ei — éz) -ns = 0. 


61 


5 Comparison of grain boundary models 


This results in Approach II (coupled slips) 


N N 
OUsY. - ðps\ 
5 (et ns- o2) t- (& nsa — 0, 


U U (5.8) 
((&& -&z) ns & (1- Qi ,ns)))=0, YaBetl...,N}, 
(5.9) 


with GB misorientation parameter £9. Examples for the misorientation 
function can be found in Bayerschen et al. (2016). In this work, the 
misorientation function is adopted from the free energy of the GBs ac- 
cording to Gurtin (2008), given in Eq. (5.1), as 


f(Qi;,ns)- C]; = (az dj) (nz x ns) - (ni x ns) . (5.10) 


5.3 Simulation setup: Shearing of a periodic 
laminate 


The gradient crystal plasticity framework, introduced in Chapter 4, is 
applied to an elasto-plastic laminate microstructure in order to discuss 
the differences between the two GB flow rules. For sake of simplicity, 
the exact solution of the model problem is discussed under the assump- 
tion of single-slip and monotonic loading. For the cyclic deformation 
of a two-phase laminate structure and the in-depth discussion of the 
resulting generalized kinematic hardening modulus, cf. Forest (2008); 
Forest and Guéninchault (2013); Wulfinghoff et al. (2015). The laminate 
is exposed to plane strain, while, per definition, the plastic shear strain 
in the elastic phase is set to zero y = 0. As shown in the schematic 
illustration of Fig. 5.1 the laminate consists of three periodic layers, two 
elasto-plastic layers A and B colored in light gray and an elastic layer E 
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a 
h 


Figure 5.1: Undeformed periodic laminate material and deformed unit cell. 


colored in dark gray. The GB between grain A and B is labeled S with 
the GB plane normal ns. The origin of the coordinate system is located 
at the GB between the two adjacent elasto-plastic layers. The widths 
of the layers A, B and E are denoted by ha, hg and hg, respectively. 
Each grain contains a single slip system, with the orientations of the 
slip direction vectors denoted by da and ds and the orientations of the 
slip plane normal vectors denoted by na and ng. The origin of the 
global coordinate system {e1, e») is located at the GB of the adjacent 
elasto-plastic grains A and B. 

In order to account for a misorientation between the two grains, a local 
coordinate system 16;, €2} is introduced in each grain, as illustrated in 
Figure 5.2. The orientation of the local coordinate systems is character- 
ized by the rotation angles o4, ag for each grain, respectively. Quan- 
tities given in the local coordinate system are denoted by a tilde, e.g., 
&. The angle between the local coordinate system and slip system is 
denoted by ®. In this work, the single slip system is not rotated to- 
wards the local coordinate system so that ? = 0. Consequently, the 
slip system corresponds with the local system yielding da = dg = € 
and na = ng = 6€». For the transformation from the global coordinate 
system {eı, e», e3} into the local coordinate system {€1, é2, 63) a trans- 
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>11 
n £2 dA e3 ni£2 ds 
€1 €1 
EE Le BY? 
ar €1 B 
" p : 
: A» ns 
! Is 
e »]« 


ha m à i” 


Figure 5.2: Notation for the orientation of the single slip systems. 


formation matrix Q, is defined by € = Q yx. The transformation is 
performed by a rotation about the e3-axis by the angles va and dp for 
grains A and B, respectively. Consequently, the transformation reads 


4 cos(0) —sin(v) 0 Ti 
i9 | =| sin(9) cos(9) 0 33 |» VI = (04, VB}. 
£3 0 0 1 x3 


(5.11) 


5.4 Analytical solution 


5.4.1 Constitutive assumptions 


The free energy in the bulk material within the small deformation frame- 
work is assumed to be additively decomposed into three terms 


(ee, ^ac; grad (y)) = Ve(&e) SEE Vs (grad (y)) a Wn(Yac) ; (5.12) 


64 


54 Analytical solution 


where ^4. denotes the accumulated plastic slip. For the exact solution a 
monotonous evolution of the plastic slip is assumed, i.e., Yac = y. The 
free energy density for the bulk thereby consists of an elastic energy 
contribution Ye, a defect energy contribution Yg and an isotropic hard- 
ening energy contribution Yh, which are introduced in the following. 
The elastic energy contribution is given by 


1 1 
QUe = 9e : Clee] = „ie € Ep) ` Cle = Ep] : (5.13) 
For sake of simplicity elastic isotropic material behavior is considered. 
Consequently, elastic material properties are independent of direction 
in space. In this case the stiffness matrix consists of only two indepen- 
dent entries 


C-AIGI-2yul, (5.14) 


where A and u are so called Lamé parameters. With this, the Hooke's 
law for linear elastic isotropic solid bodies reads 


c = Atr(e.)I + 2ue. . (5.15) 


For the defect energy density the phenomenological, quadratic ap- 
proach 


1 
e = 5 Kyllgrad (7) ||? (5.16) 
is chosen, where K, is a constant defect parameter. Isotropic hardening 


is characterized by the isotropic hardening energy 


1 
obn = zo + T07, (5.17) 
where Ko and 7 are introduced as material parameters. Supplementary 


to the free energy in the bulk material, an additional energy density per 
unit area contribution on the GB is introduced. For the special case of 
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single slip Eq. (5.1) results in 


1 
Oss = 3% (Caa Ysa + OBBYSB — 2CaBysarss) ; (5.18) 


with right-hand limit of plastic slip y* = ss, left-hand limit of plastic 
slip y7 = ysa and GB strength x. The slip-interaction moduli are given 
with Eq. (5.5) as 


CAA = (na x ns) a (na x ns) ; (5.19) 
Cop = (ns x ns): (ns X ns) , (5.20) 
CAB = (da d dg) (na x ns) s (ng x ns) y (5:21) 


5.4.2 Displacement field, strain and stress 
The displacement field for simple shear deformation is given as 
u = u (24, £2)€1 + u2(X1,X2)ea (5.22) 


in the global coordinate system. It is defined for each grain in the 
following boundaries 


UuA(r1,c2), TE [—ha, 0] 
u = į ug(z1, £2), 21 € [0, hp] : (5.23) 


ug(zi, x2) ; LE (hp, hp + hg] 


The displacement field in the local coordinate system reads 
UA, &a) 5 Tı € [—ha/ cos (Va), 0] 


41 € [0, hp/ cos (0g)] . (5.24) 
te (Z1, 4) , Ze [hp/ cos (Up), (hp/ cos (0g)) + hg] 


gr 
Il 
gi 
w 
— 
Rr 
Es 
Qu 
N 
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Considering the local coordinate system, it is further assumed that the 
displacement in €;-direction only depends on i? and the displacement 
in €>-direction only on i, yielding 


à = ü(£3)6 + fin (# Jeo. (5.25) 


An elastic isotropic and homogeneous material is assumed. The in- 
finitesimal strain tensor in the local coordinate system is defined by 


& = sym(grad (#)) = = (= dii 


= | — + — ] (€1 8 č2 +č28č1). (5.26 
2 \ dis an) (1 9 +2292) RU 
Further, with the assumption of a monotonous evolution of the plastic 
slip, the plastic strain tensor in the local coordinate system is defined 
with Eq. (3.9) by 


čp -4M* = 7 (č1 Q č2 + č2 Q 6i). (5.27) 


N| = 


The local Cauchy stress tensor is given by Hooke's law 


& = Č[ē — č] = E düz 


de da )asa+rzea). (5.28) 


The resolved shear stress in the local coordinate system results in 


= 4 die i) (5.29) 


di; dé 
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5.4.3 Analytical solution for plastic slip 


From Eqs. (4.13) and (4.54) it follows 


B= nr = Koy+ o, (5.30) 
Y 
ERO RENS 


Consequently, the bulk material flow rule, given in Eq. (4.57), results in 
Kgdiv (grad (y)) — Koy + (T — 7) = 0. (5.32) 


With Eq. (5.29) the inhomogeneous second-order linear ordinary differ- 
ential equation in the local coordinate system is given by 


d? dit, da 
i eq = -p ( — + 2 I . i 
3a 0% H ($ + da, 3) + To (5.33) 


The balance of linear momentum, given in Eq. (3.32), results in 
d 2. n 
p- Ei (E) 534, 
T5 


From this integration results in the displacement fields 


Ua = kua 2 + dua ’ (5.35) 
tip = kupt2 + dup, (5.36) 
Up = kugt2 + dur , (5.37) 
Uga = Ju diy kyAZi + dya, (5.38) 
Up = n di, + kvgřı + dvg, (5.39) 
ug = kygzi + dyz, (5.40) 
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where kua, kup, Kuz; dua; dup, dug and kya, kvg, kve, dva, dvg, dvg are 12 un- 
known integration constants. With these relations the plastic slip distri- 
bution within grain A 


736: | Ko 5 Ko- 
7a (41) = Cia exp ( Fea) + Coa exp (- ay 
5 8 (5.41) 


G (kua + kva) — To 
———— 
Ko 


and the plastic slip distribution within grain B 


"n Ko. Ko. 
(1) = cip exp yr + copexp | — y 
5 8 (5.42) 


" G (kup + ky) — To 
Ko j 


result from Eq. (5.33), with constants cı a, C24, cip, cop. In the following 
the boundary conditions are introduced to determine the 16 constants. 


5.4.4 Boundary conditions 


Boundary conditions for plastic slip 


The laminate consists of three periodic layers. Grain B, and due to 
periodicity also grain A, adjoin the elastic grain. At the respective GB 
between an elasto-plastic grain and the elastic grain the plastic slip must 
be equal zero. Therefore, the two boundary conditions read 


EE ha 

An (a ey —) =, (5.43) 
NEZ hg 
YB (a = ) =0. (5.44) 
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At the GB between the two elasto-plastic layers, the plastic slip is dis- 
continuous with a right-hand limit ysg and a left-hand limit ysa, i.e., 


^A (E 
^s ($ 


j= 
j= 


As (5.45) 


1=0 YS 
1 —0) = ss. (5.46) 
The four equations are sufficient to determine the four constants cıA, 


C2A, C1B and cop. 
Boundary conditions for the displacement field 


The laminate is exposed to a simple shear deformation. For this case 
four loading conditions can be determined. At x2 = 0 no displacement 


occurs, i.e., 


ura (21, £2 = 0) = 0, (5.47) 
urp(#1, 2 — 0) — 0, (5.48) 
uiE(z1 , £2 = 0) zu. (5.49) 


At x2 = 1 the constant macroscopic shear strain K is applied, i.e., 


UA(zi,29 = 1) = K, (5.50) 
wer, £2 = 1) = K, (5.51) 
uig(z1, T2 = 1) = K . (5.52) 


Additionally to the loading conditions the continuity of the displace- 
ment field at the GBs 


U2A (11 = 0, £2) = uap(z1 = 0, £2) 5 (5.53) 


uap(z1 = hp, x2) = uag(z1 = hg, x2) (5.54) 
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is taken into account. Eq. (3.33) states that the stress vector is continuous 
at GBs. This requires 


1(QjeQ,) ns] - o. (5.55) 


with ns = e1. Thus, one additional boundary condition is obtained for 
both the GB between the two elasto-plastic grains A and B and for the 
GB between the elasto-plastic grain B and the elastic grain E. 

Finally, uz is assumed to be a periodic fluctuation field. The periodicity 
requires 


0 hp hg+heg 
J PUB ipg J PUE dun N Quz 5i =0, (556) 
—ha 0 T h 


B Ox, 


if u2 is continuous at the GB. The definition of uz as fluctuation requires 


0 hg hg hg 
J usa day + / up dax «f u dzı = 0. (5.57) 
0 h 


—ha B 


With Eqs. (5.47)-(5.57) 12 boundary conditions are obtained that are 
sufficient to solve for the 12 integration constants of the displacement 
fields wa, up and ug. 


5.5 Comparison of grain boundary flow condi- 
tion approaches 


5.5.1 Material parameters 


The material parameters of the aluminium-type bulk material are cho- 
sen according to Bayerschen et al. (2015). Isotropic elastic material be- 
havior is homogeneous within the laminate material with = 27GPa. 
The plastic behavior in both elasto-plastic layers is also chosen to be 
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homogeneous with 7o = 10MPa, Ko = 1075MPa and K, = 84uN. The 
GB strength is chosen according to Özdemir and Yalcinkaya (2014) as 
k = 2100N/m. The remaining material parameter of GB flow condition 
approach II is chosen as £o = 100N/m. 


5.5.2 Limit case of vanishing misorientation 


In order to discuss the differences between the two GB flow condition 
approaches for the limit case of vanishing misorientation, this section 
focuses on the influence of the grain width of the elasto-plastic grains 
on the solution. In the following investigation, the width of the left 
elasto-plastic grain is chosen to be ha = 0.5um, while the width of the 
right elasto-plastic grain hg is varied between hg = ha and hg = ha/8. 
The width of the elastic layer is set to hg = 0.5um, however, with y 
being zero, the elastic layer is only shown in the figures showing the dis- 
placement field. The misorientation between the adjacent elasto-plastic 
grains A and B is assumed to vanish. Therefore, the grain structure is 
a single-crystal problem and the GB between the elasto-plastic grains A 
and B should not be observable. 

For a macroscopic shear of K — 0.01, Fig. 5.3a shows the exact solution 
of the plastic slip 5(24) in shear direction over the laminate for the GB 
flow condition approach I. In case of identical grain widths ha = hp, the 
distribution of the plastic slip is continuous across the GB. As soon as 
the grain widths of the elasto-plastic grains are not equal, i.e., hg < ha, 
the plastic slip is discontinuous with ysa > yss. It can be seen that with 
decreasing grain width hp, the discontinuity of the plastic slip at the 
GB increases. The single crystal consistency check is satisfied only for 
ha = hg, where a continuous distribution of the plastic slip is observed. 
For the limit case of a vanishing misorientation, the discontinuous dis- 
tribution of the plastic slip for the case of hg < ha is unexpected. 
According to GB flow condition approach I, given in Eq. (5.7), the gra- 
dient of the slip depends only on the GB energy ys. For the limit case 
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0.025 
0.020 - Fr, En 
— 0.015 hB = ha /2 “ii 
8 
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(a) Solution I 


1 
—0.5 —0.25 0 0.25 0.5 


xı in um xı iN um 


(b) Solution I.1 (c) Solution I.2 
Figure 5.3: Comparison of GB flow condition approach I (top) and approach II (bottom) 


for the limit case of vanishing misorientation with a macroscopic shear of K = 0.01. The 
exact solutions for the plastic slip 4(21) are shown for a varying grain width hp. 
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of vanishing misorientation the GB flow condition reads 


£g: ns +r (Ysa — yss) ^0, (5.58) 


Er ns t &(ysA — yss) = 0. (5.59) 


A solution is obtained for £4 : ns = £s: ns, i.e., the slip gradients in 
normal direction have to coincide at the GB. This can be observed in 
Fig. 5.3a for all hg < ha. 

Now considering approach II to derive the GB flow condition, Fig. 5.3b 
and Fig. 5.3c show the exact solutions of the plastic slip y(x1) in shear di- 
rection plotted over the laminate. In contrast to approach I, approach II 
results in two solutions for the distribution of the plastic slip. Solu- 
tion II.1 in Fig. 5.3b is similar to the solution obtained from approach I 
in Fig. 5.3a. In contrast to solution II.1, the distribution of the plastic 
slip in solution IL2 is not only continuous across the GB for the case of 
ha = hg, but also for all hg < ha. The single crystal consistency check 
is therefore satisfied. 

Approach II for the GB flow condition reads for the limit case of vanish- 
ing misorientation 


(£g ns + «(Ysa — YsB))YsB — (£A: ns +s (ysa — vss))*sa = 0, 
(5.60) 


£g: ns — Ea ns — 0. 
(5.61) 


Eq. (5.60) is quadratic in the plastic slips of the GB. Hence, solving 
for the two unknowns ysa and ysg results in two solutions obtained 
by approach II. The additionally introduced flow condition, given in 
Eq. (5.9), reads that the jump of the slip gradient at the GB correlates 
to a certain misorientation function, which is zero for the limit case of 
vanishing misorientation. Consequently, the slip gradients in normal 
direction of both sides of the GB coincide in both solutions II.1 and I.2, 
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see Eq. (5.61). The discontinuous solution of the plastic slip in Fig. 5.3b, 
however, fails the single crystal consistency check. 

In order to investigate differences in solutions II.1 and II.2 the displace- 
ment field u»(z1) for both solutions II.1 and II.2 is shown in Figure 5.4. 
With the boundary conditions Eqs. (5.53) and (5.54) the distribution of 
u»(z1) for solution II.1 as well as for solution II.2 is continuous at the 
GBs. The derivatives of the fluctuation field, i.e., Ou2/02, are shown 
in Fig. 5.5. The displacement field obtained by GB flow condition ap- 
proach I shows a kink at the GB in case of hg « hA. The same kink 
is observed for solution II.1 of approach II. In solution IL2, however, 
the distribution of w2(x1) is continuously differentiable at the GB. The 
differentiability of the displacement field is included in the Hadamard 
compatibility condition, which is investigated in the following. 


Hadamard compatibility condition check 
The Hadamard compatibility condition is given by 
{F}Ps — 0, (5.62) 


with GB projector Ps = (I — ns & ns), cf., e.g., Gurtin et al. (2010). For 
the simple shear of the laminate material the displacement field is given 
as u = Kage, + u2(x1)e2. From this the displacement gradient can be 
specified by 

dus 


H — K(ei ® e2) + ——(ea e) e1) A (5.63) 
dai 


Therefore, the jump of the deformation gradient yields with the assump- 
tion of small deformations 


IF] = (= u 4) (e2 G ei). (5.64) 
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hg — ha 
hp =ha/2 ------ 
hg — ha/3 ---- 
hg =ha/4 === 


(b) Solution II.1 (c) Solution I.2 


Figure 5.4: Comparison of GB flow condition approach I (top) and approach II (bottom) 
for the limit case of vanishing misorientation with a macroscopic shear of K — 0.01. The 
exact solutions of the fluctuation fields u» (x1) are shown for a varying grain width hg. 
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hp — ha 

hg =ha/2 -=-= 
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0.010 | 
— 0.005 
S 
= 0 
—0.005 
0.010 
—0.5 
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Figure 5.5: Comparison of GB flow condition approach I (top) and approach II (bottom) 
for the limit case of vanishing misorientation with a macroscopic shear of K — 0.01. 
The exact solutions of the derivatives of the fluctuation fields 9u2/0x, are shown for 
a varying grain width hp. 
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For the considered three-phase periodic laminate, the normal vector of 
the GB is given by ns = ei. The Hadamard condition reads 


d d 
[F]Ps = ( de = za) (e & e1) (ea & ea + es & es) 


| 


Cras a) ((e1 : €2) (e2 B ea) + (e1 - ex) (es & ex) 
0. 


(5.65) 


Considering the case of no grain misorientation, the Hadamard com- 
patibility condition is therefore trivially fulfilled. Hence, even though 
the displacement is continuously differentiable at the GB only for so- 
lution II.2, the Hadamard compatibility condition is not violated for 
solution IL1. 


5.5.3 Misaligned grains 


In the following, the influence of a varying grain misorientation is in- 
vestigated for GB flow condition approach I and approach II. In order 
to investigate the influence of the grain misorientation, a variation of 
the misorientation angle Ja = {0°, 10°, 20°, 30°} is chosen. The three 
layers of the laminate are equally sized with ha = hg = hg = 0.5um. 
Figure 5.6a shows the exact solution of the plastic slip plotted over the 
laminate structure with the use of GB flow condition approach I. The 
laminate is exposed to a macroscopic shear strain of K = 0.01. For the 
equally sized grains without misorientation the single crystal solution 
is obtained. For all solutions including a grain misorientation, the distri- 
bution of the plastic slip is discontinuous at the GB between the elasto- 
plastic grains. It can be seen that with increasing grain misorientation 
angle, the discontinuity of the plastic slip at the GB increases. 

With approach II for the GB flow condition, again two solutions for 
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Figure 5.6: Comparison of GB flow condition approach I (top) and approach II (bottom) 
for misaligned grains with a macroscopic shear of K = 0.01. The exact solutions of the 
plastic slip 7(%1) are shown for a varying misorientation angle dx. 
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for misaligned grains with a macroscopic shear of K = 0.01. The effective stress-strain 


Figure 5.7: Comparison of GB flow condition approach I (top) and approach H (bottom) 
curves are shown for a varying misorientation angle A. 
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the plastic slip are obtained. The exact solutions II.1 and II.2 of the 
plastic slip are shown in Fig. 5.6b and Fig. 5.6c. Both solutions show 
a discontinuous distribution of the plastic slip for misaligned grains, 
more distinct than obtained by approach I. In contrast to the investiga- 
tion of a varying grain size, neither of the two solutions is identical to 
the solution obtained by GB flow condition approach I. Equivalent to 
approach I, however, it can be observed that with an increasing misori- 
entation angle, the discontinuity of the plastic slip at the GB increases. 
Additionally, based on the amount of the grain misorientation, the gra- 
dient of plastic slip increases. This can be related to the amount of 
dislocations, which are piled-up at the GB. 
In order to investigate the differences of the solutions obtained by GB 
flow condition approach I and approach II, the laminate is exposed 
to an increasing macroscopic shear strain of K € {0,...,0.01}. This 
allows the investigation of the GB states during loading. In Fig. 5.7, the 
effective stress-strain curves are shown for both GB flow condition ap- 
proaches. While for approach I, shown in Fig. 5.7a, only two sectors of 
the stress-strain curve are observable, three GB states can be identified 
for approach II, which is shown in Figs. 5.7b and 5.7c: 

* Elastic: No yield condition is satisfied, therefore, no plastic deforma- 
tion occurs. The plastic slip is zero. 

* Microhard: The yield condition is satisfied only in the bulk material 
but not in the GB. Dislocations pile-up. 

* Microcontrolled: The yield condition is satisfied in the bulk material 
as well as in the GB. GB plasticity occurs, which is controlled by a GB 
flow rule. 

For GB flow condition approach II the two solutions, IL1 in Fig. 5.7b 

and IL2 in Fig. 5.7c, show an equivalent stress evolution. The slope of 

the microhard GB stress evolution is higher for a higher misorientation 
angle of the adjacent slip system. The most dominant influence of the 
misorientation angle, however, is the transition point of the microcon- 
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trolled GB state. The resistance of the GB against dislocation transmis- 
sion strongly depends on the misorientation function, cf. Eq. (5.9). For 
Ja = 0°, i.e., limiting case of coinciding slip systems, dislocations can 
cross the GB without hindrance. Consequently, the GB is microfree after 
the initial elastic state. The amount of dislocations piled-up at the GB 
before the critical stress for GB transmission is reached is higher for a 
higher misorientation angle. 

To illustrate the different GB states and the differences between ap- 
proach I and approach II for the GB flow condition, the distribution of 
the plastic slip along the z-axis is shown in Fig. 5.8 for the four loading 
states labeled in Fig. 5.7. Exemplary, the slip distributions are shown 
for Ya = 30°. For approach I shown in Fig. 5.8a plastic slip evolves 
at the GB as soon as the initial critical resolved shear stress occurs. 
Approach II again shows two solutions shown in Fig. 5.8b and Fig. 5.8c. 
It is important to note that only for the microcontrolled GB state, two 
solutions are obtained. Both solutions are identical for the microhard 
GB state, where the GB still acts as an insurmountable obstacle for 
plastic slip. After a sufficient amount of dislocations piled-up at the 
GB, the GB flow condition is fulfilled and GB slip is permitted. For this 
microcontrolled GB state the discontinuity of the plastic slip evolves 
and two distinct solutions are obtained. 


5.5.4 Investigation of size effect 


In order to investigate the size effect of the periodic bicrystal the 
laminate width h = ha = hg = hg is varied between h = 0.25um and 
h = 10um. Traditionally, the experimental dependence of the strength 
of a polycrystal is described by the inverse square-root of grain size 
via the Hall-Patch relation, cf. Hall (1951). Li et al. (2016), however, 
investigate various datasets from literature and conclude, that their 
data strongly supports a simple inverse or logarithmic dependency. 
Arguments are made against the Hall-Patch approach, whereas the 
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Figure 5.8: Exact solutions for the evolution of the plastic slip 5(21) for GB flow condition 
approach I (top) and approach II (bottom) with grain misorientation Ya = 30°. 
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simple inverse and logarithmic dependency are supported by un- 
derlying dislocation theories. Following Li et al. (2016) the elastic 
strain Ee = 7/G, i.e., stress normalized by the shear modulus, and 
normalized width d = h/ao, i.e., laminate width normalized by a lattice 
constant ao, are relevant parameters. The lattice constant is chosen 
to ag — 0.361nm. Consequently, the following three size-dependence 
equations are investigated 


ks In (d) 
d 


_ = 1 = = 2 _ M 
£e] = E01 + — Em = E0 +— , £el3 = E03 + , (5.66) 


Vd’ d 
where £01, £02, £03 and k1, k2,k3 are the corresponding fitting parame- 
ters. The parameter fit is performed with the nonlinear least-squares 
fitting procedure in MATLAB. 

Fig. 5.9 shows the fitting approaches of Eq. (5.66) for the dependence 
of elastic strain on the normalized laminate width for K = 0.01 based 
on GB flow condition approach I (left) and GB flow condition approach 
II (right). The corresponding root-mean-square errors are listed in Ta- 
ble 5.1. All three fitting approaches show a good agreement with the 
analytical results, however, the Hall-Patch approach of 1, shows a 
higher deviation from the analytical results, especially for bigger grain 
sizes. The simple inverse size-dependence equations, i.e., E12, predicts 
the analytical results best. 

The analytical results obtained with GB flow condition approach I and 
approach II show a similar size-dependency. For small grain sizes, how- 
ever, the individual volume-averaged stress results differ. For larger 
sizes the behavior of the bulk material dominates the mechanical re- 
sponse of the periodic three-phase laminate. For small sizes, however, 
the GB flow condition approach stronger influences the mechanical re- 
sponse. The comparison of £,; for equal d shows higher strains obtained 
by GB flow condition approach II. 
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Figure 5.9: Dependence of the elastic strain £j = 7/G on the normalized laminate width 
d= h/ag for K = 0.01 based on GB flow condition approach I (left) and GB flow 
condition approach II (right). The analytical results are compared to the three fitting 
approaches of Eq. (5.66), exemplary shown for 94 = 0° (red) and Ya = 30° (green). 


Table 5.1: Root-mean-square errors for size dependency of the elastic strain £j = 7/G 
according to the three approaches of Eq. (5.66) with GB flow condition approach I and II. 


ApproachI 
EN T. s 


Approach II 
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5.6 Summary and conclusion 


With the gradient crystal plasticity model the dislocation transmission 
through the GBs is investigated, which is of special interest during the 
plastic deformation of metal materials. For the modeling of the GBs 
two approaches are introduced. The common approach to model the 
interaction of dislocations with GBs as introduced by Gurtin (2008) is 
based on the assumption of independent plastic slip rates at the GB. It 
is shown that with the commonly made assumption of independent 
plastic slip rates at the GB the consistency check for coinciding slip 
systems is failed. In order to avoid this shortcoming, a second approach 
to derive the GB condition is introduced. In this new approach, an addi- 
tional flow condition is applied, where the discontinuity of the plastic 
slip gradient at the GB correlates to a misorientation function. Applying 
the new approach to the limit case of vanishing misorientation results in 
two different solutions. While one solution is an artificial, non-physical 
solution with strong similarity to the solution obtained by the common 
approach, the second solution shows a single-crystal response with a 
continuous and symmetric parabolic distribution of the plastic slip even 
under variation of the grain size. 

Additionally to the consistency check for the limit case of vanishing 
misorientation, the analytical solutions are compared for misaligned 
grains. While under assumption of independent plastic slip rates at 
the GB only a microcontrolled GB state is observed, an additional mi- 
crohard GB state is identified for the new GB flow condition approach. 
This microhard state results in a pile-up of dislocations at the GB until 
the introduced flow condition is satisfied. After the flow condition 
is satisfied in the bulk as well as in the GB, transmission across the 
GB occurs and a discontinuous distribution of the plastic slip evolves. 
Consequently, the two stages of GB behavior which are observed in Sun 
et al. (1998), i.e., buildup of dislocations near the GB and a vanishing 
peak at higher strain levels, can be modeled. Both, the critical pile- 
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up stress for dislocation transmission as well as the amount of slip 
discontinuity across the GB depend on the GB flow condition, i.e., on 
the misorientation angle between the two grains. 

For both GB flow condition approaches the effect of the laminate width 
on the mechanical response of the periodic three-phase laminate is in- 
vestigated. The analytical results show a good agreement with a sim- 
ple inverse or logarithmic approach, rather than the traditional inverse 
square-root Hall-Patch approach. This result is supported by Li et al. 
(2016) where arguments are made against the Hall-Patch approach and 
the simple inverse and logarithmic approach are supported by underly- 
ing dislocation theories. 
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Chapter 6 


Accumulated plastic slip 
framework 


6.1 Introduction and motivation 


The following chapter is based on Erdle and Bóhlke (2017). A gradient 
crystal plasticity theory is introduced in a thermo-mechanical frame- 
work based on the principle of virtual power. Hereby, the micromor- 
phic approach, as developed in Forest (2009), is used. In the previous 
sections for each slip system one gradient-stress is introduced in order 
to account for long range dislocation interactions. Consequently, each 
slip system requires one additional DOF in the FE implementation for 
the nonlocal continuum model. Considering all plastic slips massively 
increases the computational cost compared to a classical continuum me- 
chanical description. To reduce the DOF in crystal simulations, Wulfin- 
ghoff and Bóhlke (2012) introduce an accumulated plastic slip gradient 
plasticity model. For details on the computational benefit regarding 
this approach it is referred to Wulfinghoff et al. (2013). Whereas the 
model preserves the single crystal slip kinematics, it is restricted to one 
gradient-stress, associated with the gradient of the accumulated plastic 
slip, in order to account for long range dislocation interactions. While 
the gradient of the plastic slip on an individual slip system has a direct 
physical interpretation in terms of the density of geometrically neces- 
sary dislocations, the gradient of the accumulated quantity can only 
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be interpreted as an approximated measure of the overall dislocation 
pile-up structure. 

Experimental studies, e.g., Lee et al. (1990), have shown discontinu- 
ous plastic slips on individual slip systems across GBs. Whereas the 
accumulated plastic slip is assumed to be continuous in Wulfinghoff 
et al. (2013), the presented model accounts for discontinuities of the 
accumulated plastic slip at GBs. The GB flow rule is evaluated at sharp 
interfaces using discontinuous trial functions in the ABAQUS FE imple- 
mentation. This approach is based on the theory of the extended FE 
method, as presented in Fries and Belytschko (2010). 


6.2 Model approach 


6.2.1 Basic assumptions 


Taking the plastic slip on each slip system into account, one additional 
DOF for each slip system is required. The computational cost of the 
resulting strongly coupled system of generally non-linear equations is 
massively increased, compared to a classical continuum mechanical de- 
scription. In order to reduce the DOF in crystal simulations, a simplified 
gradient crystal plasticity theory based on an accumulated plastic slip 


2N . 2N 
^fac = Y falar - MAI (6.1) 
a=1 a=1 


is introduced in Wulfinghoff et al. (2013). Hereby, for N slip systems 2N 
directional dependent plastic slips Aa are introduced. The slip parame- 
ters are non-decreasing by definition, i.e., Àa > 0. The plastic slip of a 
slip system is given by the difference between its two directional depen- 
dent plastic slips, i.e., ya = Aa — Aa+n. Although formally a physical 
simplification may be associated with a framework based on an accu- 
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mulated plastic slip, the numerical efficiency is increased notably. The 
simplification by use of the accumulated plastic slip was tested in vari- 
ous simulations, e.g., Wulfinghoff and Bóhlke (2012); Wulfinghoff et al. 
(2013), where it showed a good agreement with size effects observed in 
experiments. 

Following the work of Forest (2009), the micromorphic variable ^4. as- 
sociated with yac is introduced as an additional DOF in order to account 
for microstructural effects. The concept of introducing a micromorphic 
medium to account for micro-structural effects is based on the early 
work of Eringen and Suhubi (1964); Mindlin (1964). The additional 
DOE, its gradient and the deviation of the micromorphic variable to its 
associated counterpart contribute to the power of internal forces. For an 
overview of the micromorphic concept and the general procedure of the 
enhancement of a continuum to incorporate generalized micromorphic 
gradient effects it is referred to Forest (2009); Forest et al. (2011). 


6.2.2 Principle of virtual power 


Based on Eq. (4.24) the following contributions of the bulk material and 
the GBs are assumed to define the microscopic system. The bulk mate- 
rial contributions are given by the Cauchy stress ø that expends power 
conjugate to the velocity gradient L, by a scalar micromorphic stress 7 
conjugate to Jac and by a vector valued micromorphic stress € conjugate 
to the rate of the micromorphic slip gradient grad (5,.). In addition to 


internal power in the bulk material, GB micromorphic tractions II^ are 
assumed to expend power conjugate to the left-hand limit +, and the 
right hand limit 4. of the micromorphic accumulated plastic slip rate 
at the GBs. Consequently, the power of the internal forces reads 


Pint = f (a: L+ Tac + €- grad (Jac)) dv 
p (6.2) 


+f (II + II.) da. 
S 
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An alternative form of Eq. (6.2) is obtained with the use of the jump 
ata GB {4ack = 37, — Jac in combination with the mean value at a GB 
(jac) = (Fx + Vac)/2. This description of the discontinuity provides 
an easier interpretation and comparability to existing theories with a 
continuous solution field, i.e., {act > 0, (ac) — Jac: By using the 
jump and the mean value of the micromorphic accumulated plastic slip 
rate, the power of the internal forces reads 


Pint = / (e Lc T jac + £ ` grad (Žac) ) dv 
y 


+ f[ (n 0.) da, 


(6.3) 


where the work conjugate quantities of the jump and the mean value 
(I* -II), — Il, 2 II* - II- (6.4) 


are introduced. The power of external forces is constructed based on 
Eq. (4.25). It is assumed that the macroscopic system is defined by a 
body force density ob and a macroscopic surface traction t, both con- 
jugate to à, and by a micromorphic surface traction = conjugate to “ac. 
The corresponding parts of the external boundary are denoted by 3Y, 
and OYz, respectively. Consequently, the external power has the form 


Pow = | eb-üao+ f Ludos | SA... da. (6.5) 
y ay OVs 


The principle of virtual power is constructed in a similar manner to 
Section 4.3.1. With the application of divergence theorem the power 
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balance reads 
f (m — div (€)) b4ac dv — / (ob + div (o)) - ù dv 
y y 
= (II; 15 Yac} + Im (öYac) m 1£654.] j ns) da (6.6) 


+f (en =E) Sheeda + f (en — t) -óu da — 0, 
OYz Ovi 


where the virtual solution fields {6%, 6/4. } can be varied arbitrarily. The 
following local field equations, Neumann conditions and GB condition 
are obtained 


r — div (€) =0, Va € Y, (6.7) 

gb + div (o) — 0, Vrcy, (6.8) 

Th, {Yack + Ha sc) — [£344] ns = 0, Va € S, (6.9) 
£-n—E-0, Va € Vz, (6.10) 

on-t=0, Va € 0y,. (6.11) 


6.2.3 Helmholtz free energy 


The mass specific Helmholtz free energy of the bulk material ~ is split 
into four contributions 


y= We( Fe) T Wh (Nac) Ar V; (Grad (Fac)) sr Vx (Vac = Yac) , (6.12) 


i.e., elastic energy v., isotropic hardening energy wv», gradient defect 
energy v, and penalty energy Yy, which is of numerical nature. 


A Hooke-type elastic contribution is considered 


1 
doe = zEe CIEe], (6.13) 
where the reference mass density oo is introduced. 
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Isotropic hardening effects of the bulk material are modeled through 
the phenomenological hardening energy 


1 
Qon = z Oac + (Teo x T0 )^fac 
(Qo mE E ei) 


Too — TO 


iso (Too — 70)” exp (- 
cf., e.g., Palm (1951). Hereby, the initial and asymptotic values of the 
critical resolved yield stress 7) and Tə as well as the initial and asymp- 
totic values of hardening modulus Oo and ©, are introduced as mate- 
rial parameters. For monotonic loading processes, it can be shown that 
a direct connection exists between the Voce-type bulk material hard- 
ening energy and the dislocation-density-based Taylor formulation of 
hardening, cf. Bayerschen (2016). 
Additionally, a gradient defect energy contribution introduces an inter- 
nal length scale. Examples of defect energy contributions are given in 
Reddy (2011). Based on the small deformation framework of Wulfin- 
ghoff et al. (2013), a quadratic energy contribution based on the gradient 
of the micromorphic accumulated plastic slip is used 


d T 
09g = 5 KsllGrad (ac) I? (6.15) 


with the defect parameter K, to account for energy of geometrically 
necessary dislocations. A non-quadratic choice of a power-law defect 
energy based on the gradient of the accumulated plastic slip is investi- 
gated in Bayerschen and Bóhlke (2016) for a single-crystal gradient plas- 
ticity theory. Note, that in contrary to a defect energy which depends on 
the sum over gradients of plastic slip, discussed, e.g., in Reddy (2011), 
Eq. (6.15) depends on the gradient of the sum over plastic slips. 


Finally, based on the work of Forest (2009), a numerical penalty energy 
1 N 2 
QoUx = 3 Hx ac — Yac) (6.16) 
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is introduced. The penalty energy penalizes deviations of the micro- 
morphic accumulated plastic slip Ya. from the associated counterpart 
^ac by a large penalty factor H, and consequently ensures the coupling 
of both variables. 

An additional energy contribution is assumed to arise from the GB re- 
sistance to plastic flow. The GB energy is assumed to depend on the 
jump and the mean value of the micromorphic accumulated plastic slip. 
In the following, the special case of a linear GB energy is considered 


esos = Kj{¥ack + Bol Yau) y (6.17) 


with reference GB mass density oso and material parameters K; and 
Km. The phenomenological approach of Eq. (6.17) accounts for a discon- 
tinuous distribution of the micromorphic accumulated plastic slip. This 
is an extension of the linear GB energy based on the accumulated plastic 
slip used in Wulfinghoff et al. (2013). As a consequence of the linear 
dependency of the GB energy on the jump and the mean value of the 
micromorphic accumulated plastic slip, GB hardening is neglected. The 
hardening extension of the GB with a quadratic GB energy contribution 
is investigated in Bayerschen et al. (2015). Since the parameters K; 
and Km do not depend on the misorientation of slip systems at GBs, 
orientational effects are neglected. This is a major difference to the GB 
energy of Gurtin (2008) investigated in Section 5.2. 


6.2.4 Clausius-Duhem inequality and flow rules 


The dissipation inequality for the bulk material and the GBs is given 
with Eq. (4.35) by 


/ oo) dV. < Pexto(Vo), Í esos dA < Pexto(So) , (6.18) 
Yo 


So 
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where Pexto(Vo) denotes the external power expended on the bulk ma- 
terial in the reference configuration and Pexto(Vo) denotes the external 
power expended on the GBs in the reference configuration. By a trans- 
formation of Eq. (6.3) the internal energy contribution with respect to 
the reference placement is obtained 


Pinto = n (S : E + T0 fac + &o . Grad (Nac) ) dV 
Vo 


(6.19) 
+f (Tjo {Vac} + Imo Yan) ) dA . 
0 
The stresses in the reference placement read 
S=det(F)F'oF', m=det(F)r, £,-—det(F)F !£, von 
Ilo = det(F)Il, IL,;o—det(F)l,. 
Bulk material 
With Eqs. (6.13)-(6.16) if follows 
j OWe r ape r On apy : 
yY = -E+ : F, + ( + Yac 
OE OF, Mac Mac (6.21) 


OV, z Oz 


toe Oad 


- Grad (Fac) ' 


Inserting the derivation of the Helmholtz free energy into the Clausius 
Duhem inequality of the bulk material, given in Eq. (6.18), results in 


= O00%e ; = Joov>~ M 
D= (s JE ) E+ (n O5, "yas 


Joow : 
+ (& -= i - Grad (Jac) (6.22) 


ooe ; Joon , OpoUx Y. 
— * == —— B Pa a A, > 
OF, P ( Mac n Mac ee 9s 
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where the power balance and Eq. (6.19) are used. The assumption of 
purely energetic stresses S, 7 and £, results in the potential relations 


_ Oe. OE. _ pag -T 
S= 5p, gp Fr EF, (6.23) 
oop y 
To = TM = HN Yoc) (6.24) 
O 09s "m " 
= CCS. ac * 2 
Šo Grad (4c) K,Grad (Yac) (6.25) 


Furthermore, the reduced dissipation inequality of the bulk material 


Dass RO, F,- Ut + Rx) Asc > 0 (6.26) 
Mac Mac 


is obtained. The following notations are introduced for the derivative 
of the isotropic hardening energy density and the penalty energy by the 
accumulated plastic slip 


Mac Mac 


Consequently, Eq. (6.26) can be written as 
Drea = Me i Ly m (8 +p) ^fac 2 0, (6.28) 


where the Mandel stress M, — FIF.S. is introduced. Based on the 
Mandel stress the resolved shear stress of slip system a is given as 
Ta = Me- Ma. Eq. (6.28) motivates the bulk material flow rule 


DENT u (6.29) 


TD 


where {-} denotes the Macaulay bracket with property {-} = max (-,0). 
Hereby, the reference shear rate Ào, the initial critical shear stress 7o, the 
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drag stress Tp and the slip rate sensitivity p are introduced as material 
parameters. 


Grain boundaries 


With Eqs. (6.17), (6.18) and (6.19) the GB dissipation per unit area reads 


a Oosots Y ‚: 
Dee (n. = z— [Sach + (n. - AS (Fac) 20. (6.30) 


The energetic GB micromorphic tractions are given as 


e _ O0soVs è Josovs 


{Fach (ac) 


With the GB condition, given in Eq. (6.9), and under the assumption of 


dissipation free GBs, the GB dissipation inequality results in 


0 = (K,(Grad (Jac) * nso — Kj) [jac] 


E n (6.32) 

+ (K,{Grad (ac) “TSO — Kn) (se) : 

The rate independent GB yield criteria are specified by 
(pj = K,(Grad (Yac)) “NnS — K; (6.33) 
Ym = K,{Grad (Fac)} - nso — Km. (6.34) 


From this the loading conditions for Eq. (6.32) can be formulated in the 
Kuhn-Tucker form with 


Pj ack = 0, Yj < 0, {Fact = 0, (6.35) 
Vis ac) = 0, mS 0, (Fae) = 0. (6.36) 
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6.3 Finite element implementation 


6.3.1 Enrichment of shape functions 


In the following the enrichment of FEs for discontinuous solution fields 
across GBs is discussed. In order to model a discontinuous displace- 
ment field of crack simulations, an enrichment of the classical FE ap- 
proximation is presented in Moés et al. (1999). A generalization of 
the model for arbitrary discontinuities, called extended FE method, is 
introduced in Belytschko et al. (2001). An alternative approach to the 
enrichment of FEs would be the cohesive interface approach, based 
on the pioneer works of Barenblatt (1959). This approach is widely 
used for a continuum mechanical theories, e.g., in Gurtin and Anand 
(2008) for dislocation-based inelastic processes by slip and separation 
at inter-crystalline GBs. The extended FE method used in this work, 
however, has the advantage that the FE mesh does not need to conform 
to the polycrystal topology, cf., e.g., Simone et al. (2006). This allows 
the modeling of various crystalline structures or evolving interfaces 
without remeshing. 

The micromorphic accumulated plastic slip field is approximated by the 
standard FE ¥-(@) = N (x) and the additional nodal unknowns to 
account for discontinuities. Therefore, the general structure of shape 
function enrichment reads 


N? (£) 
= 
Fac(®) = I, N? (s); X NP (a) (dw) -ó(m))à, —— (637) 
icly icly 
standard FE enrichment 


where Iy is the set of all nodes and ¢(x) denotes the enrichment func- 
tion, cf., Fries and Belytschko (2010). The enrichment function is evalu- 
ated at the integration points as well as at the nodes of the FE x;, Vi € Iy. 
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Due to the choice of the enrichment term, the standard FE approxi- 
mation is retained for FEs which do not contain a GB. Based on the 
choice of the enrichment function, the solution field is either decoupled 
completely (strong discontinuity), or the solution is continuous while al- 
lowing for kinks across interfaces (weak discontinuity). Typical choices 
for the enrichment function for both weak discontinuities and strong 
discontinuities are given in Fries and Belytschko (2010). In this work the 
strong discontinuities of the accumulated plastic slip field are modeled 
with the sign of the level-set function 


’ wer) ns «0; 
ó(r) 240,  (x—zms):ns —0, (6.38) 


; (z— cs): ns » O0. 


In general, the level-set function is capable to model dynamic surfaces, 
cf., Osher and Fedkiw (2001). In the work at hand, however, GBs are 
assumed as fixed material singular surfaces. 


6.3.2 Discretization of grain boundaries 


In the following a mesh-independent GB discretization is introduced. 
The construction of grain structures is done by means of the grain- 
indices of the FE nodes, as illustrated in Fig. 6.1. Consequently, image 
data, e.g., Electron Backscattered Diffraction scans can be used as a 
starting point for the simulation. If the grain-indices of a FE vary, a 
GB is contained which is assigned to integration points depending on 
predefined cases. Because the structure within a FE is prescribed, the 
area of the GB inside the FE can be partitioned and allocated to the 
integration points. Additionally, it is ensured that within one FE, all 
values of the enrichment function, given in Eq. (6.38), are present. Con- 
sequently, the solution field discontinuities are modeled in agreement 
with Eq. (6.37). An additional DOF is introduced, which represents the 
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Node of grain 1 


Node of grain 2 


* 


Integration point of grain 1 


* 


Integration point of grain 2 


* 


Integration point of grain boundary 


Figure 6.1: Exemplary distribution of nodes and integration points for a three- 
dimensional eight-node finite element containing a GB. 


grain index of each node. In doing so, the local orientation of the GB is 
obtained by determining the gradient of this indices-field. The gradient 
of the indices-field is nonzero if and only ifa GB is contained inside the 
FE. In that case, the orientation of the GB can be obtained by a mapping 
of the additional nodal index DOF with the shape function gradients 
onto the integration points. The normalization of this gradient results 
in the direction of the strongest variation of the grain indices which can 
be defined as the normal vector of the GB. Consequently, it is possible 
to implement curved or polyhedral surfaces, although the mesh size 
influences the quality of the surface approximation. As a result of the 
information about the local orientation of the GB and the discretization 
of the surfaces with the use of the integration points, a direct allocation 
and distribution of surface energy is enabled. The GB flow rule, given 
in Eq. (6.32), is evaluated accurately at the integration points. 
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6.3.3 Newton scheme 


Local Newton scheme 


The theory at hand is implemented using the user element interface of 
ABAQUS CAE. Equivalent to the standard FE procedure with elasto- 
plastic user materials, the local update of stresses is performed until 
defined residuals are reached. In this work the update of the Cauchy 
stress is covered by updating the elastic deformation gradient. With 
Eqs. (3.14) and (3.15) the first residual equation for the local Newton 
scheme reads 


2N 
rı =F,- FR — At (er. -F.), inti) . (6.39) 


a=1 


The subscript n denotes the current time step. The subscript n + 1 for 
quantities of the subsequent time step is dropped. For the general- 
ized relative stress p, Equations (6.7) and (6.24) imply the additional 


residuum equation 
r 2N 
ra= Mac + qp — Yee — At Kal. (6.40) 
x a=1 


Eq. (6.39) does not fulfill the necessary condition det(F.) = det(F') re- 
sulting from an isochoric plastic deformation, automatically. After an 
updated F% fulfilling Eqs. (6.39) and (6.40) is found, the determinant 
correction is performed 


off det(F) V? 
P= Fr: (a) (6.41) 


For consistency with the algorithmic implementation an array-matrix 
notation is used in the following, where 2nd-order tensors are repre- 
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sented as arrays, e.g., rı > rı. In order to determine the algorithmic 
tangent, the total derivatives of the residuum equations, given in 


Eqs. (6.39) and (6.40), are set to vanish, i.e., 


| Ory Or, 
Ore Org 
= dE —— dö+ — 
at 


In matrix-vector notation this reads 


Ora Sri 

OF, dp |[9E- 
a ir ; 
ar, Op) \@ 


Or, Ory, 4. 
OL [dL] w Mac Fac , (6.42) 
Ore Ore " 
aL Et ace m 
of ii Or 1 
ƏL WIE. 
(6.44) 
m. ‚DB a 
OL Fac vee 


As a result, the algorithmic tangent for 


given by 
dE. GE, Or, 
dL ac OF. 
dp dp Ore 
dL ac OF e 


Global Newton scheme 


the local Newton scheme is 


A 


The computation of the field variables u = N" ü and 44, = N” TN d ö 
on the global level within a FE framework, requires the use of aNewton 
scheme with associated algorithmic tangent. The weak form of the 
balance of linear momentum, given in Eq. (6.8), is obtained with a trial 
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function dw" = N" w by an integration over an arbitrary volume Y 


0 = 7 div (c) - dw" dv, (6.46) 
v 


pe n B"'cdv- | N"'tda. (6.47) 
v Yi 

Hereby, body forces are neglected and the matrix of shape function gra- 

dients B" is introduced. The weak form of the micromorphic field equa- 

tion, given in Eq. (6.7), is obtained with the trial function dw? = N’ ù, 

GB condition and GB flow rule as 


0 = f (B — p — Kadiv (grad (Ya.))) dw? dv, (6.48) 
y 


= I NIT (8 — d) dv «f B''K, (B74 + 2°) dv 
2 ¥ (6.49) 


_ NZa | Ky da, 
OYa Ss 


act 


where Sact denotes the plastically active part of the GBs. The enrich- 
ment of the accumulated plastic slip by Eq. (6.37) results in an addi- 
tional equation 


j= l (8 — B— K,div (grad (5,:))) dw dv, (6.50) 


= | xS @-pavt | BP i (me ii) d 

» " d (6.51) 

- N? Edo « f Kjda, 
Ove S; 


act 


with the trial function óu? = N° ®. The residuum equations, given in 
Eqs. (6.47), (6.49) and (6.51), are solved for {û, 7, ô } by a Newton scheme 
within ABAQUS CAE. Details on the global algorithmic tangent, which 
is required for this procedure are given in Appendix A.2. 
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0.25 um 


25 um 


Figure 6.2: Illustration of boundary conditions of laminate microstructure (top) and 
representative bicrystal (bottom) with periodic boundary conditions. 


6.4 Simulation setup 


FE simulations are performed for a laminate microstructure, illustrated 
in Fig. 6.2 (top). By the use of periodic boundary conditions, a rep- 
resentative bicrystal with 25um x 0.25um x 0.25um in dimension, as 
illustrated in Fig. 6.2 (bottom), is investigated. 

At first, the bicrystal is loaded displacement controlled along the z- 
axis up to an overall shear strain of 0.05. In order to account for an 
approximated stationary slip case, an adaptive time stepping algorithm 
is used during the simulation time of 10°s. An active set search is 
performed on the GBs in order to control the convergence of the global 
Newton scheme. The global Newton scheme is considered to be con- 
verged, if the residual tolerance is sufficiently small, and if the set of 
plastically active GB integration points has not been changed during an 
iteration. For details on the active set search on the GBs it is referred to 
Wulfinghoff et al. (2013). 

A FCC crystal structure is considered. The stiffness tensor is charac- 
terized by the elastic constants C1111 = 168GPa, C1122 = 121GPa and 
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Table 6.1: Model parameters for simple shear simulations of periodic bicrystal. 


0MPa 330MPa 330MPa 


| 


C1212 = 75GPa. The slip directions and normal vectors of the slip sys- 
tems are listed in Appendix A.1. The mismatch of slip systems at the GB 
is realized by a misorientation of adjacent grains. One grain is oriented 
in a standard (100)-orientation with respect to the x-axis of the Carte- 
sian coordinate system. The adjacent grain is rotated by various angles 
between o; = 0° and yy = 90? along the x-axis. The material parame- 
ters of the bulk material and GBs are varied in order to study param- 
eter influences on the mechanical response of the crystalline structure. 
If not further specified the material parameters are chosen based on 
Wulfinghoff et al. (2013), where the model parameters for simulations 
of copper microwire tensile tests were fitted to the experimental data of 
Yang et al. (2012). The material parameters of this work are listed in Ta- 
ble 6.1. The FE mesh is structured and composed of user-defined, cube- 
shaped hexahedrons with 27 integration points. For the DOF of the 
displacement field, quadratic shape functions are used. All additional 
DOF which represent plastic strain fields or are used for numerical 
reasons (e.g., indices field), are interpolated with linear shape functions. 
A convergence study is performed for the simple shear loading of the 
laminate microstructure with mesh discretizations between 10? DOF 
and 10? DOF. The mismatch of adjacent grains is induced by a grain 
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Figure 6.3: Convergence study for a varying mesh discretization of the crystalline 
structure. A good compromise of the computational time and the accuracy of the results 
is given for a FE mesh with about 6 x 10* DOF. 


rotation angle of v, = 45°. In Fig. 6.3 the influence of a varying mesh on 
the final average true stress o,. and on the final average accumulated 
plastic slip yac is shown. A good compromise of the computational 
time and the accuracy of the results is given for a FE mesh with about 
6 x 10* DOF. Consequently, this discretization is used for the following 
simulations. 


6.5 Finite element results 


6.5.1 Influence of grain misorientation 


In Fig. 6.4 (left) the resulting stress-strain curves are shown for various 
grain misorientation angles. The three domains of the shear stress evo- 
lution, i.e., elastic, microhard and microcontrolled are recognizable. The 
first domain, until the initial critical shear stress To is reached, the shear 
stress only depends on the shear modulus which is almost identical for 
each grain structure. The critical yield stress at which plastic deforma- 
tion of the bulk materialis initialized, however, varies for different crys- 
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tal orientations. As soon as the bulk material deforms plastically and 
while the GB still acts as an insurmountable obstacle for plastic slip, the 
hardening of the bulk material is a combination of linear isotropic hard- 
ening and the hardening induced by dislocation pile-ups. This domain 
is strongly influenced by the misorientation angle. Even though the 
misorientation of neighboring grains is not contained in the GB energy 
in a natural way, the misorientation angle strongly influences the shear 
stress which allows for dislocation transition across the GB and conse- 
quently initializes the third domain. In contrary to Section 5.2, however, 
the GB does not vanish for a grain misorientation angle yy = 0°. The 
misorientation of neighboring grains and the orientation of the GBs 
is not contained in the GB energy in a natural way. Admittedly, no 
GB would exist between adjacent grains, nonetheless, a resistance of a 
virtual boundary to the passing of dislocations is considered. A sep- 
arate consideration of orientational effects is required for a consistent 
theoretical framework. As soon as the GB flow rule is fulfilled and 
the loading condition is satisfied, GB slip is permitted. Due to linear 
dependency of the GB energy and the herewith resulting neglection 
of GB hardening, the subsequent stress evolution is barely influenced 
by the misorientation angle of adjacent grains. Additional dislocations 
arising inside the bulk material are allowed to glide across GBs without 
additional hindrance. 

Fig. 6.4 (right) shows the corresponding final distributions of the micro- 
morphic accumulated plastic slip over the x-axis of the grain structure 
for various grain misorientation angles yy. Due to the symmetry of 
the periodic bicrystal only one half with a GB located at x = 6.25um 
is shown. For equally oriented grains, i.e., Yx = 0°, a continuous and 
symmetric distribution of the micromorphic accumulated plastic slip is 
obtained. The virtual GB in-between both grains acts as an obstacle 
until a critical amount of dislocations piled-up and the critical state 
fé}- ns = Kj, (£) ns = K is reached. This critical state is not depen- 
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Figure 6.4: Influence of grain misorientation angle i»; on the mechanical response (left) 
and corresponding distribution of the final micromorphic accumulated plastic slip Yac 
over the length of the periodic bicrystal (right). 


dent on the misorientation angle between the two grains, consequently, 
the jump and mean value of the slip gradient in normal direction of 
the GB are the same for various v; values. The discontinuity of the 
micromorphic accumulated plastic slip, however, strongly depends on 


the misorientation angle. 


6.5.2 Size effect for equally oriented grains 


For the following investigations, tensile test simulations are performed 
for the periodic microstructure, illustrated in Fig. 6.2. The bicrystal is 
loaded strain controlled up to an overall strain of 0.1 along the z-axis. 
An adaptive time stepping algorithm is used during the total simulation 
time of 100s. Both grains are considered to be oriented equally in a 
standard (100)-orientation with respect to the x-axis of the Cartesian 
coordinate system. In order to investigate the influence of model param- 
eters while neglecting orientational effects, this special case of p, = 0? is 
considered. As a consequence of equally oriented grains, a comparison 
of the numerical results to existing investigations with a GB energy 
independent of the discontinuity of the accumulated plastic slip at the 
GB, e.g., Wulfinghoff et al. (2013), is possible. In the following simula- 
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0.25um lm M m "E NN | n=5 


Figure 6.5: Illustration of simulated grain sizes. The grain size is modified by a variation 
of the number of laminates n inside the 25um periodic structure. 


tions, the GB energy parameter associated with the discontinuity of the 
accumulated plastic slip is neglected, i.e., K; = 0Nm~!. The remaining 
model parameters are chosen as listed in Table 6.1. As illustrated in 
Fig. 6.5, the size effect is investigated by the variation of the laminate 
width inside the periodic structure. 

Fig. 6.6 shows the influence of the number of laminates n inside the 
25um periodic structure on the mechanical response. For a decreas- 
ing grain size, the model responds with an increased strength. Conse- 
quently, size-effects occurring on the microscale are reproduced. The 
elastic range of the bulk material is independent of the investigated 
grain size. The hardening of the bulk material induced by dislocation 
pile-ups, however, is strongly influenced by the grain size. After GB 
slip is permitted, the stress evolution is identical for the different grain 
sizes. Due to the identical orientation of adjacent grains, additional 
dislocations arising inside the bulk material are allowed to glide across 
GBs without hindrance. 
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Figure 6.6: Influence of the number of equally sized grains within a crystal structure of 
identical width on the mechanical response. As illustrated in Fig. 6.5, a higher n correlates 
to a smaller grain size which results in a higher stress response of the periodic bicrystal 
under tensile load. 


6.5.3 Parameter influences for equally oriented grains 


In the following, the influence of the GB energy parameter Km and the 
defect energy parameter K, on the mechanical response of the peri- 
odic bicrystal under tensile load is discussed. In Fig. 6.7 the GB en- 
ergy parameter Km is varied in a range between Km = 0Nm-! and 
Km = 200Nm-1, while all other parameters are kept constant. In Fig. 6.7 
(left), the stress-strain curves are shown. For Km = 0, no resistance of 
the GB to the passing of dislocations exists. Consequently, no disloca- 
tions pile-up at the GBs. As shown in Fig. 6.7 (right), the corresponding 
distribution of the final micromorphic accumulated plastic slip is ho- 
mogeneous over the x-axis of the grain structure. Upon an increase 
of Km, the second domain of the stress-strain curve in which plastic 
deformation occurs in the bulk material and the GBs still act as insur- 
mountable obstacles for plastic slip extends, leading to a stiffer material 
response. The effect of the GB energy parameter on the amount of 
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dislocations which pile-up at the GBs before the GB flow rule is fulfilled, 
is clearly evident in Fig. 6.7 (right). As soon as GB slip is permitted, the 
stress evolution only depends on the multiplication and annihilation 
effects of statistically stored dislocations. These effects are modeled by 
the isotropic bulk material hardening energy which is identical for the 
different GB energy parameters. Consequently, the evolution of the 
stress-strain curve in the third domain is independent of the choice 
of Km- 

In Fig. 6.8, the parameter of the defect energy K, is varied between 
K, = 10 4N and K, = 1N, while the GB energy and all other param- 
eters are kept constant. Note that the choice of K, = ON is not valid 
due to the evaluation of the GB yield condition. The defect energy 
parameter K, only shows an influence on the second domain of the 
stress-strain curve. After the GB flow rule is fulfilled, the stress-strain 
curves coincide. Upon an increase of K, the pile-up of dislocations at 
GBs is penalized with an increased energy contribution. Consequently, 
a higher value of K, results in a more homogeneous distribution of 
the final micromorphic accumulated plastic slip ^4. over the x-axis of 
the grain structure. For K, — 1N, the GBs still act as insurmountable 
obstacles for plastic slip at a true strain of 0.01. 

Finally, the influence of the isotropic hardening energy parameter Oo 
is illustrated in Fig. 6.9. By increasing the isotropic hardening energy, 
the slope of the stress-strain curve in the third domain is increased. At 
the same time, the slope in the second domain is slightly increased, as 
well. Physically, these dependencies are explained by the multiplication 
and annihilation effects of statistically stored dislocations which occur 
during the whole plastic range. The effect of dislocation pile-ups at 
the GBs, however, dominate in the second domain. The corresponding 
distributions of the micromorphic accumulated plastic slips are shifted 
along the y-axis to lower slip values, while the influence of the GB is 
hardly affected. 
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Figure 6.7: Influence of the GB energy parameter Km on the mechanical response (left) 
and corresponding distribution of the final micromorphic accumulated plastic slip ^ac 
over the length of the periodic bicrystal (right). 
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Figure 6.8: Influence of the defect energy parameter Kg on the mechanical response (left) 
and corresponding distribution of the final micromorphic accumulated plastic slip ^ac 
over the length of the periodic bicrystal (right). 
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Figure 6.9: Influence of the isotropic hardening energy parameter Oo on the mechanical 
response (left) and corresponding distribution of the final micromorphic accumulated 


plastic slip Yac over the length of the periodic bicrystal (right). 
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6.6 Summary and conclusion 


The presented gradient crystal plasticity, finite strain theory of Erdle 
and Bóhlke (2017) allows the modeling of dislocation phenomena on 
the microscale of the material. Instead of considering slip gradients on 
individual slip systems, a gradient-stress associated with the gradient 
of the micromorphic accumulated plastic slip accounts for long range 
dislocation interactions. The use of the discontinuous trial function in 
the FE implementation allows the evaluation of a GB flow rule in due 
consideration of the misorientation of adjacent grains. In doing so, an 
energetic contribution of GBs controls the pile-ups of dislocations, with- 
out the evaluation of a GB dislocation density tensor. In contrary to the 
GB energy considered in Chapter 5, the misorientation of neighboring 
grains and the orientation of the GBs is not contained in the GB energy 
in a natural way. A separate consideration of orientational effects is 
required for a consistent theoretical framework. 

Different types of energy contributions of the micromorphic, discontin- 
uous field variable and its gradient are modeled. The FE implementa- 
tion is discussed as well as the influence of the energy parameters on the 
mechanical response and corresponding distributions of the plastic slip 
over periodic laminate structures. The presented approach provides a 
promising numerically efficient framework for large scale FE simula- 
tions which includes the grain size dependence of the yield strength. 
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Transport mechanisms and grain 
boundaries 


7.1 Introduction and motivation 


Recent studies by Hochrainer et al. (2014) introduce the production of 
the total dislocation density dependent on line curvature as well as the 
ability to account for dislocation transport within a higher-dimensional 
continuum dislocation dynamics framework. Based on this framework 
a thermodynamic formulation is constructed in Hochrainer (2016). By 
taking the kinematics of dislocations by continuous field variables into 
account, a physically based theory of crystal plasticity is obtained. Sin- 
gle crystal continuum dislocation dynamic simulations show a good 
agreement with the experimental studies of micro-bending tests and the 
discrete dislocation dynamic simulations of Motz et al. (2008). Schulz 
et al. (2017) show by a comparison with discrete dislocation analyses 
that the continuum description can capture dislocation interaction pro- 
cesses in a heterogeneous materials. A first three-dimensional numeri- 
cal implementation with a discontinuous Galerkin scheme is presented 
in Sandfeld et al. (2015). Wulfinghoff and Bóhlke (2015) formulate a 
three-dimensional FE implementation of the dislocation continuum the- 
ory within a gradient crystal plasticity framework. On the basis of the 
gradient of an accumulated plastic slip, the model is still coupled to 24 
partial differential evolution equations of the total dislocation densities 
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in a FCC crystal. In the following chapter an approach is presented in 
order to regain a physically motivated framework while preserving the 
computational benefit of accumulated field variables. 


7.2 Dislocation field evolution 


7.2.1 Equation of motion 


In the following dislocations are considered on a scale where they can 
be identified individually. Nevertheless, a coarse-grained description 
is obtained again in a phenomenological sense (cf., e.g., Monavari and 
Zaiser (2018)). Dislocation climb and dislocation cross-slip is neglected. 
The motion of a dislocation on a slip system o is considered as a propa- 
gating front moving normal to itself with velocity vy, i.e., the velocity is 
orthogonal to the line direction ta of the dislocation at each point. This 
implies that the line direction is on the corresponding slip plane o, i.e., 
ta ` Nna = 0. Consequently, the drag force fa of a dislocation (per unit 
length) on a slip system a, is related with its velocity va by a viscous 
drag coefficient Ba via the relation 


fa -Bava = —Bava(na x ta), (7.1) 


cf. El-Azab (2000). The evolution of a dislocation line is caused by 
two major effects. On the one hand, the application of external stresses 
and stresses generated by other crystal defects results in a force acting 
on a dislocation line. On the other hand, however, the curvature of 
a dislocation line can also be considered as a contribution to the local 
internal stress of the dislocation, cf. Mughrabi (2001). Consequently, the 
velocity of the dislocation line may be a function of the self-interaction 
caused by the line curvature and of an underlying flow generated by 
external stresses and by other crystal defects. 
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In Osher and Sethian (1988) it is shown that equations of motion of prop- 
agating fronts with such velocity dependencies resemble Hamilton- 
Jacobi equations. Let I(t) be a family of curves, generated by the 
evolution of a reference curve along an in-plane vector field with 
velocity vr. The evolution equation of propagating fronts constructed 
in Osher and Sethian (1988), the so called level-set equation, is based on 
the level-set function ór(z, t). The shape and position of a propagating 
front is uniquely determined as the zero-set of the corresponding 
level-set function: T(t) = {x | ör(z,t) = 0). A level-set framework for 
the movement of a curve by curvature is constructed by Ambrosio and 
Soner (1996), where the mean curvature is introduced as 


kr = -div (grad (¢r) ||grad (¢r) ||") . (7.2) 


Note that grad (ór) is ambiguous for a given zero-set of the level-set 

function, however, the mean curvature kp is uniquely defined. The 

level-set equation, constructed by Osher and Sethian (1988), results as 
dpr 


p X vr ||grad (er) || = 0. (7.3) 


For a detailed discussion of the level-set method it is referred to the 
review article of Osher and Fedkiw (2001). 


7.2.2 Local form of the dislocation transport 


After a single dislocation has moved through a volume of the crystal 
on a slip system o, a distortion of the crystalline lattice remains. The 
absolute value of the Burgers vector b defines the distance of slip caused 
by the dislocation. The angle y describes the angle between the line 
direction and Burgers vector of the dislocation. By the introduction of 
this angle ọ as an additional dimension, dislocations can be described 
in a higher-dimensional configuration space. The total dislocation den- 
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a) c) * 


Figure 7.1: a) Discrete dislocation network on slip system o, b) Discrete distortion of the 
crystal lattice along line A-A with dislocation velocity va, c) Continuous plastic slip along 
line A-A with velocity of plastic slip field v3. 


sity pt (æ,t) is expressed in terms of a orientation distribution func- 
tion po (a, t, v) for dislocations 


f )of* 
pi (x, t) = x ], Qa (2,0, 9) do, (7.4) 
cf. Hochrainer (2013). With the assumption of an isotropic dislocation 
velocity v, (z,t,q), ie. vo(z,t) = Valæ, t, p), Hochrainer et al. (2014) 
construct an evolution of the total dislocation density based on a statis- 
tical dislocation density measure. The statistical measure is restricted 
to dislocation networks where all dislocations form circular loops or 
straight lines, cf. Hochrainer et al. (2007). 

In this work, instead of the construction of an evolution equation based 
on the higher-dimensional configuration space, continuous plastic slip 
fields G4 representing the distortion of the crystal lattice caused by a 
network of discrete dislocation lines are considered. For an illustration 
of this transfer, see Fig. 7.1. As illustrated in Xiang (2009) the density 
distribution of discrete dislocation lines on each slip plane is identified 
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by contour curves of a dislocation density potential function after local 
homogenization. Zhu and Xiang (2015) introduce an additional poten- 
tial function in order to account for the distribution of the slip planes 
in the averaged dislocation continuum. In this work, however, it is 
sufficient to describe the plastic slip on each slip plane and identify 
the density distribution of dislocation curves in the averaged disloca- 
tion continuum. This approach results in a connection between the 
net length per volume dislocation density on each slip system over a 
representative volume and the in-plane gradients of plastic slip. Ashby 
(1970) introduces the density of GND p$ based on the absolute value of 
the in-plane gradient of plastic slip as 


||Grad? (Ca) || = pad, (7.5) 


where Grad? (6,) = (I — ng & n4) Grad (Ga) denotes the in-plane gra- 
dient of plastic slip. In Monneau and Patrizi (2012) Hamilton-Jacobi 
equations which describe dislocations by the Peierls-Nabarro model, 
cf. Nabarro (1997), are homogenized in order to construct non-local 
equations for densities of dislocations on the macroscopic scale. This 
study motivates the application of the level-set approach of Osher and 
Sethian (1988) to the evolution of the micromorphic plastic slip 


Wa 


P 
Da = ys, Grad? (ee 


à Pe I Grad? (Ca) || . (7.6) 


Hereby v7, denotes the velocity of the continuous plastic slip field (see 
also Fig. 7.1). The combination of Eqs. (7.6) and (7.5) results in 


a Gp * 


"Ot PaVa (7.7) 
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By taking the time derivative on both sides of Eq. (7.5) it follows 


Ope 


P 
75 = -Div (pv) — pGvzDiv eed . (78) 


Grad" (Ca) || 


The first part of the right-hand side of Eq. (7.8) accounts for dislocation 
transport, while the second part accounts for curvature induced pro- 
duction by the curvature 


ee Grad? (C4) 
u an: 09 


cf., e.g., Xiang (2009). The basic equation of the local dislocation trans- 
port is obtained 


Ops 


a =D (oSv) + pgviks. (7.10) 


7.3 Coarse-grained dislocation transport 


7.3.1 Details on the coarse graining operator 


Let V; be a suitable representative volume element whose centroid is 
located at the origin. Points in V. are denoted by Z. The coarse graining 
operator 


V(X,t) = (V(Z,t).:— z g(X —Z(t))U(Z,t)dV (711) 


with the weighting function g(a) is introduced. The weighting function 
is normalized, smooth and vanishes smoothly at points close to 0Q¢, i.e., 


1 
— X)dV 21, lim X)=0. 7.12 
y f, 200 im, aO) (712) 
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The weighting function is assumed to be isotropic. For details on the 
coarse graining operator, it is referred to section 2.6 in Nemat-Nasser 
and Hori (1993). In the following, dependencies of quantities are 
dropped due to better readability. By introducing the coarse-grained 
densities of dislocations p, = (p$ )., the time derivative results in 


Opa 1 Og G Ops 


Ue 5 y, oe +9 8t dV. (7.13) 

It follows 9 8 oZ da AZ 
c MON en (7.14) 

ot O(X—Z) at OX Ot 

Consequenty, the first part of Eq. (7.13) results in 
Og G u Og Oz g 
cl, a zh- gr^» V 

° (7:15) 


1 Oz 1 Oz 
— Deel gee Ze Go? 1A. 
Ve J, i (ws a) ud Velen OF 


The flux integral over the area of the representative volume element O9. 
vanishes. It follows that the order of applying the coarse graining oper- 
ator and differential operators Div (-) and Grad (-) are interchangeable, 
i.e., (Div (V)), = Div ((V).) and (Grad (V)). = Grad ((V).). 


7.3.2 Semi-phenomenological dislocation transport 


Eq. (7.10) is coarse-grained by means of a representative volume Qe. 
This volume is assumed to be time independent. Coarse grained quan- 
tities are denoted by V = (W).. The coarse-grained form of Eq. (7.10) 
results in P 

Te = —Div ((pSut) ) + (oSv; k;), . (7.16) 
The velocity vz is assumed to be homogeneous inside Qe, i.e. v; = Da, 
cf. Hochrainer et al. (2014); Monavari and Zaiser (2018). Consequently, 
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Eq. (7.16) results in 


“Pe = —Div (para) + (oh, ; (7.17) 
with the coarse-grained velocity 


=~ - UL) D 
Va = —— = — Grad? a 
Wate — E (Grad? (Ca) 


as (7.18) 
The coarse-grained in-plane gradient of plastic slip on slip system a is 
given by 

(Grad? (6,)) 2 PaGrad ((¢a)c) . (7.19) 


Applying the coarse graining operator to Eq. (7.7), the Orowan equation 
is obtained 


c3 = ai = (pn bvi) , = PabPa . (7.20) 
Micromorphic plastic slips Ya = (Ca). carry the microstructure effect, 
i.e. the targeted gradient effects. As discussed in Forest (2009), a close 
connection exists between the micromorphic quantities and the corre- 
sponding macro counterparts, i.e. Yq. 

In this work an approach is presented in order to adopt the above phys- 
ically motivated framework to preserve the computational benefit of 
accumulated field variables. The numerical implementation by Wulf- 
inghoff and Bóhlke (2015) requires 24 dislocation field variables. As 
opposed to this high computational cost, in this work only two addi- 
tional field variables, i.e., the accumulated density of dislocations and 
the accumulated mean curvature, are introduced in order to account 
for the dislocation transport of Eq. (7.10). Under the assumption of 
proportional loading with respect to each slip system, Eq. (7.10) can be 
expressed in terms of accumulated plastic slips on each slip system Yaca 
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and the corresponding in-plane gradient as 
j 1 x G,*1* 
rad Caca) ) + (punks). (7.21) 


This simplification allows the construction of a model based on the 
accumulated plastic slip similar to Chapter 6. Note that this assumption 
does not imply monotonic loading processes, implicitly, but monotonic 
evolutions of plastic slip on each slip system. As a consequence, if the 
slip direction of a slip system changes its sign, for example due to a 
cyclic loading process, no annihilation of dislocations occurs. 


In order to formulate the model based on the accumulated plastic slip, it 
is necessary to formulate an approach connecting the in-plane gradient 
of the slip on slip system a with the in-plane gradient of the accu- 
mulated plastic slip. As an approach the gradient of the accumulated 
plastic slip is prorated on the individual slip systems by 


O' laco, v Po Mac y = 
P, (acra (Yaza) = z. gi Grad (Fac) = 0, (7.22) 


Pac 


where the accumulated density of dislocations fac = >>, Pa is intro- 
duced and Pa = (I — Nna & n4) denotes the projector onto the slip sys- 
tem a. The approach preserves the slip kinematics of in-plane glide of 
dislocations on slip systems. It is ensured that the direction of in-plane 
gradient of slip is preserved. Note, that the sum of the gradients of 
plastic slip on the individual slip systems is only equal to the sum of the 
gradient of the accumulated plastic slip if 5 >, Bia ‘aca Grad (Vgc) 2 0 
is valid. By using this approach it is possible to construct a numer- 
ically efficient micromechanical framework based on accumulated 
field variables which accounts for the transport of dislocations on a 
coarse-grained scale. Whereas in the model of Wulfinghoff and Bóhlke 
(2015) 24 partial differential evolution equations of the total dislocation 
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densities in a FCC crystal are introduced, only one additional DOF is 
needed in the following. 

By combining Eqs. (7.7), (7.21) and (7.22), and using $^, Pa = SI, the 
evolution of the density of dislocations on slip system a results in 


N 


Grad (fae) ) + DO G8), 023) 


a=1 


Pac Op: Mac 
ðt = Div (IE, at p ^ 


In the following the coarse-grained curvature induced line length pro- 
duction is neglected inside the representative volume element V. 


N 
6 ee (7.24) 


Small velocities are assumed which results in the identity of the partial 
time derivative and the material time derivative (cf. Eq. (3.2)). With the 
introduction of the accumulated dislocation velocity 


TM E (7.25) 
Pacb 


the semi-phenomenological dislocation transport equation is obtained 
" ; 1 , 
pac — Div (rad D . (7.26) 


7.33 Principle of virtual power 


Two adjustments are made in comparison to Eq. (6.2). The vector val- 
ued micromorphic stresses £, is introduced conjugate to the rate of the 
in-plane micromorphic accumulated plastic slip gradient grad? (4ac). In 
addition, the micromorphic accumulated plastic slip rate is assumed be 
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continuous at GBs. Consequently, the power balance is is postulated by 
N . 

I (= — x div a) ac dv — I (ob + div (o)) - ù dv 
v peri v 
N 

= I Patat ns | Hac da 

L(r- Bones) 

+f (3 asm - 8) Fina f (en — t) -óüda — 0, 

av= av, 


(7.27) 


where the virtual solution fields {$ü, 05.) can be varied arbitrarily. The 
field equations, GB condition and Neumann conditions read 


div (c) = 0, Vrcy, (7.28) 
N 
7 — M div(P4£,) — 0, Va € y, (7.29) 
asl 
N 
E- [Y Pata} ns - 0, Vx € S, (7.30) 
a=] 
on—t=0, Vx € OV, (7.31) 
N = 
MI hE =U, Va € Vz. (7.32) 
a=1 


7.3.4 Helmholtz free energy 


The mass specific Helmholtz free energy of the bulk material ~ is split 
into three contributions: 


N 
Y = pel Fe) + M Vn (Pac, Grad? (Yaca)) +x (YacrYac), (7.33) 


a=l 
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i.e., the classic, elastic energy we, the hardening energy wv; and the 
penalty energy yy, which is of numerical nature. 


Following Section 6.2.3, a Hooke-type elastic contribution is considered 


Loe = SB. $ CIE.] . (7.34) 


Zaiser (2015) identifies three contributions to the energy of a dislocation 
system. A non-local contribution based on the density of GND wv, a self- 
energy contribution 7; and a correlation contribution pe. With these 
contributions the hardening energy is given as 


Vn — Vg T Ws F We . (7.35) 


In the work of Hochrainer (2016) small densities of GND are assumed. 
Consequently, the non-local contribution based on the density of GND 
is neglected, however, this term is bound to stay, cf. Zaiser (2015). 

For complex dislocation structures explicit expressions of the line ten- 
sion of a dislocation, the self-energy per unit length, as constructed in 
Hirth et al. (1966), are not applicable. For such cases Hirth and Lothe 
(1982) give an approximation for the line tension of a dislocation by 
T ~ ub?. Based on this relation the non-local energy is introduced in 
Ohno et al. (2008) based on the density of GND. A gradient theory for 
single-crystal plasticity using this approach is derived from the princi- 
ple of virtual power by Gurtin et al. (2010); Gurtin and Ohno (2011). 
For monotonic loading processes this approach is consistent with the 
work of Ortiz and Repetto (1999); Hurtado and Ortiz (2013), where the 
non-local energy is formulated based on the absolute values of in-plane 
gradients of plastic slips in a linear form. In this work a non-local 
energy contribution is introduced in a quadratic form 


N 
Kg j 
Botts = pb? ($ V |IGrad? (ac) 9 (7.36) 


a=1 
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with the energy parameter K,. This form is based on the previous 
gradient plasticity theories dependent on the accumulated plastic slip, 
e.g., Erdle and Bóhlke (2017); Wulfinghoff et al. (2013); Bayerschen et al. 
(2015). Note that despite the formulation of Eq. (7.36) based on the ac- 
cumulated micromorphic plastic slip, the non-local energy contribution 
is introduced by means of in-plane slip gradients in order to account 
for the in-plane glide of dislocations on slip systems. The energy con- 
tribution is proposed using the Lagrangian gradient. For a discussion 
regarding alternative formulations using the Eulerian gradient and the 
gradient with respect to the intermediate configuration it is referred to 
Ling et al. (2018). 

Based on the work of Groma et al. (2007), where the dislocation energy 
is obtained by a statistical theory of dislocations, and motivated by 
Zaiser (2015); Hochrainer (2016), the self-energy 


doth. = ub? (-io« In (&)) (7.37) 
PO 
and the correlation energy 
K N 
00 Uc = pb? (z P 5 [Grad? (Fac) 9 , (7.38) 
ac mdi 


with the constants of order one K, and K. are introduced. As discussed 
by Hochrainer (2016) the reference dislocation density po is not a scaling 
parameter, but a well defined physical quantity, namely, the saturation 
value for the density of admissible dislocations. A comparison of the 
energy contribution to the dislocation energy derived by Wilkens (1969) 
results in pọ = b~?. 
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Finally, following Section 6.2.3, a numerical penalty contribution based 
on Forest (2009) is introduced in the form 


Lu vs 
009 = 5 Hy jac — tac)" ; (7.39) 


with the penalty modulus H,. By use of a large penalty modulus, the 
coupling of the accumulated plastic slip to its micromorphic counter- 
part is ensured. 

For resistance of GBs to plastic flow the GB energy approach of Wulf- 
inghoff et al. (2013) is used. Hereby, the assumption of a continuous 
accumulated plastic slip field across the GBs is made. The GB energy 
is assumed to depend on the micromorphic accumulated plastic slip a 


linear form 
0soUs = Km^jac , (7.40) 


with a GB energy parameter Km. Microhard and microfree GBs are rep- 
resented by the two limit cases Km — oo and Km = 0 respectively. Asa 
consequence of the linear dependency of the GB energy, a hardening of 
the GBs is neglected once GB slip is activated. 


7.3.5 Clausius-Duhem inequality and flow rules 


The Clausis-Duhem inequality is constructed in a similar manner to 
Section 6.2.4 with the internal energy contribution with respect to the 
reference placement 


N 


Pinto = / (s E oae + X. Eao > Grad? Gn) dV 
Vo (7.41) 


a=l 


+ 1 Moac dA. 
So 
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Bulk material 


With Eqs. (7.34)- (7.39) the derivation of the Helmholtz free energy of 
the bulk material is given by 


= Ode ; Oe Oy + Oy : 
ab = DE E+ IF Fy + Lr Pac + ac Yac 
T N (7.42) 
Xv g P (xg 
um Yac + dGra d? 5 a) Grad (Fac) * 


Inserting into the Clausius Duhem inequality of the bulk material, 
given in Eq. (6.18), results in 


E O obe B Ooo» M 
D= (s -5E ) E+ (n — Du. Nac 
N 
Joon $i 
2 — ————— ] - 7.4 
+ e Grad? = Grad (Fac) ( 3) 


ooe ; Oeo p 1 ’ Ooo» . 
..900Ve | E > 
OF, Fy - “Bi D Vac „Grad (^ac) Dou ^ac > 0, 


where the dislocation transport equation, given in Eq. (7.26), is used. 
The divergence theorem is applied to the last part of the volume inte- 
gral, resulting in 


deoe ~ Joo i 
p= (sag) +>: (eo sna) on Go 


8 JooWn " Joo» Y e 
+ (n + ENT. Grad (er) $ Grad (^ac) — a Yac 
O00Ve p  Doovy . 
= .F a > 
ap ^g em 
(7.44) 
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Additionally, the following Neumann condition is obtained 


doon 
Pac 


5.5 Grad (Yo) =0,  VeeóV. (7.45) 


The assumption of purely energetic stresses S, mo and ao results in the 
potential relations 


_ Jooe OE. piñ -T 
Ee on d C[E,]F;' , 
O 0o Un " " 
oo ar (Iu Pac ) “Grad (Tie) = P, (7.46) 
O oon Kep 2 
Su o: PIS K, ub P(x). 
Eao 7 3Grad? (ac) (= Das Grad i) 


Furthermore, the reduced dissipation inequality of the bulk material 
Drea = Me: Ly m p ^fac 20 (7.47) 
is obtained. The herewith motivated bulk material flow rule 


oe p 
NER [em (7.48) 


TD 


strongly resembles Eq. (6.29). Due to the missing energy contribution 
dependent on the accumulated plastic slip, however, no hardening con- 
tribution of £ is included. 


Grain boundaries 


With Eq. (7.40) and the divergence theorem which resulted in Eq. (7.45) 
the GB dissipation per unit area is obtained 


z Joovn 8 Josovs \ : 
—[2 ac — ace ac 2 0. 
Ds ( o + | TA FE) Grad (X ) $ ns EI Y 0 
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The following notation is introduced for the derivative of the hardening 
energy density by coarse-grained densities of dislocations 
_ Oooh 


Bp = Din (7.50) 


In this work, dissipative effects are neglected and the GB condition, 
given in Eq. (7.30), has to be fulfilled. The rate independent GB yield 
criterion is specified by 


N 
. Kop 2 8 Pry 
y= P» ( od Kaub” + a) Grad? (u) ns — Km (7.51) 


The loading conditions read in the Kuhn-Tucker form 


("jac = 0, p< 0, Jac = 0); (7.52) 


7.4 Finite element implementation 


The enrichment of shape functions, discretization of GBs as well as 
the local newton scheme is consistent with Section 6.3. On a global 
level, however, the equations for the bulk material to compute the field 
variables for any given time increment At differ. 


* Balance of linear momentum with neclection of body forces 
J div (c) dv = 0 (7.53) 
v 


e Auxiliary equation for 8p 


f (4+ Ka? (m (BE) +1) + E graa ao) l) dv =0 
y po Pac 
(7.54) 
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* Field equation for micromorphic accumulated plastic slip 


f Hy (Nac — Yac) dv = / div (s pet + Ku?) grad cs) dv 
y y ac 


8 ? 
+f z, 18d (B5) - grad (jac) dv 
(7.55) 


* Field equation for dislocation transport 


/ Inc doc / div (rad D) du 255) 
y v Pach” 

In total the ABAQUS user-element has eight DOF: Three for the dis- 
placement field u, one for the auxiliary quantity £,, one for the micro- 
morphic accumulated plastic slip Yac, one for the accumulated density 
of dislocations Jac and an additional grain indices-field in order to im- 
plement curved or polyhedral GBs. For details on the FE discretization 
and the computation of the local and the global Newton scheme it is 
referred to Section 6.3. 


7.5 Simulation setup 


For the FE simulations a FCC crystal with N = 12 slip systems is used. 
The elastic material behavior is characterized by the three elastic con- 
stants C1111 = 108GPa, C112» = 62GPa and C1212 = 28.3GPa. As dis- 
cussed in Wulfinghoff and Bohlke (2015), these elastic constants de- 
scribe a nearly isotropic elastic material behavior with a shear modulus 
of u = 26.12GPa. Based on Wulfinghoff and Böhlke (2015) the bulk ma- 
terial flow rule is characterized with reference shear rate yo = 10~°s~!, 
drag stress Tp = 1MPa and slip rate sensitivity p = 20. For the initial 
critical resolved yield stress the Taylor approach To = aotıb/ Pin: is used, 
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cf. Taylor (1934). Hereby, following the work of Wulfinghoff and Bóh- 
Ike (2015), the initial dislocation density of each slip system is chosen 
as Jini = 3001m ? and the absolute value of the Burgers vector is set 
to b = 0.234nm. With the dimensionless parameter ao = 0.3 the initial 
critical resolved yield stress results as ro ~ 32M Pa. 

The hardening energy material constants are chosen as K, = K. = 1 
and K, = (Pinib®) This choice results in a equal contribution of Yg 
and v, for p = pis. The GB flow rule is described by the GB energy 
parameter Km which is chosen as Km = 100Nm-!. The remaining, nu- 
merical parameter, is chosen sufficiently large as Hy = 10°MPa. 
Displacement controlled FE simulations are performed for a laminate 
microstructure, illustrated in Fig. 7.2 (top). By the use of periodic 
boundary conditions it is sufficient to investigate a representative 
bicrystal with dimension lym x 1.25nm x 1.25nm, as illustrated in 
Fig. 7.2 (bottom). The bicrystal is loaded along the z-axis up to an 
overall strain of 0.01. Both grains are oriented in a standard (100)- 
orientation with respect to the z-axis of the Cartesian coordinate system. 


7.6 Finite element results 


In Fig. 7.3 (top), the stress-strain curve for the tensile test simulation is 
shown. The volume-averaged true stress cz; is plotted over the volume- 
averaged true strain ln (U) = Inv F'F. Hereby, three domains are 
recognizable. Fig. 7.3 (left) shows the distributions of the accumulated 
plastic slip over the length of the periodic bicrystal for the three stress- 
strain states at times tı, t» and t5, labeled in Fig. 7.3 (top). The corre- 
sponding distributions of the accumulated density of dislocations are 
shown in Fig. 7.3 (right). Due to the periodicity of the bicrystal only 
one half of the length is shown, with the GB located at x = 0.25um. 
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Figure 7.2: Illustration of boundary conditions of laminate microstructure (top) and 
representative bicrystal (bottom) with periodic boundary conditions. 


At first the bulk material as well as the GBs deform elastically. Due to 
the absence of plastic slip in this domain the bulk material elastic-plastic 
transition point only depends on the initial critical shear stress 79. As 
soon as plastic deformation occurs inside the bulk material, dislocations 
arise inside the bulk material, transport of dislocations occurs and dis- 
location pile-up at the GBs. The GBs are still microhard in this domain, 
meaning, they act as insurmountable obstacles for plastic slip. The 
hardening of the bulk material is a combination of isotropic hardening 
of dislocations and dislocation pile-up at the GBs. After the critical GB 
state is reached, the GB flow rule is fulfilled and the loading condition 
is satisfied, dislocation transport across the GB is permitted and, conse- 
quently, a plastic deformation of the GB occurs. This activation point of 
the GB slip is mainly influenced by the self-energy parameter K,. For 
small values of K,, the resistance of the GB to the passing of dislocations 
is reduced. The corresponding distribution of the final micromorphic 
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Figure 7.3: Mechanical response (volume-averaged) of tensile test simulation (top). The 
three vertical lines correspond to the three different GB states as indicated. Distribution of 
the micromorphic accumulated plastic slip ^c (left) and distribution of the accumulated 
density of dislocations pac (right) over the length of the periodic bicrystal. The GB is 
located at x = 0.25um 
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Figure 7.4: Evolution of the volume-averaged accumulated density of dislocations pac 
and volume-averaged accumulated density of GND JS = ||Grad? (ac) |^ !. The 
three vertical lines correspond to a microhard GB, the activation of GB slip and a 
microcontrolled GB. 


accumulated plastic slip is more homogeneous over the length of the 
grain structure. Due to the linear form of the GB energy, additional 
dislocations arising inside the bulk material are allowed to glide across 
GBs without additional hindrance. Consequently, no hardening of the 
GB occurs and the resulting force on the evolution of plastic slip is 
constant. The corresponding distribution of the micromorphic accumu- 
lated plastic slip is shifted to higher slip values for increasing strain, see 
Fig. 7.3 (left). Dislocations which are piled-up at the GBs annihilate as 
soon as GB slip is permitted, see Fig. 7.3 (right). This decrease of the av- 
erage accumulated density of dislocations, shown in Fig. 7.4, results in 
a softening in the microcontrolled GB state, as shown in Fig. 7.3 (top). In 
contrary to phenomenological gradient plasticity frameworks the dislo- 
cation transport equation assures, without consideration of production 
terms, that the averaged dislocation density is constant in the micro- 
hard GB regime. As soon as GB slip is activated, piled-up dislocation 
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annihilate at the GB, which results in a decline of number of dislocations 
within the simulated grain structure. Additionally, the evolution of the 
volume-averaged accumulated density of GND 7€ = ||Grad? (^4) ||b 7! 
is shown. The results show that the portion of GND densities on the 
total density of dislocations is small, however, can not be neglected. 


Additionally to the tensile test simulation of the laminate microstruc- 
ture a simple shear simulation is performed for a fiber-matrix composite 
material as shown in Fig. 7.5. Periodic boundary conditions are applied 
to a representative segment of the composite material, which is loaded 
along the x-axis up to an overall displacement of Auo = 10nm. The com- 
posite material consists of elastic fibers of diameter 650nm embedded 
in an elasto-plastic matrix material in a periodic manner. For simplicity, 
the fibers are assumed to have the same elastic properties as the matrix 
material. In the matrix material dislocations pile-up at the GBs until the 
GB yield condition is fulfilled. Even after the activation of GB slip, how- 
ever, the boundaries still act as insurmountable obstacles for plastic slip. 
Due to the elastic behavior of the fibers plastic slip occurs only in the 
matrix material. Consequently, the pile-up structure of the dislocations 
at the GBs, i.e., the amount of GND 59 = ||Grad (Jac) ||», enlarges in 
the matrix material. Inside the fibers the density of dislocations does 
not evolve and hence has a homogeneous value of Jini. The small band 
inside the fibers with dislocation density values between Jini and the 
piled-up dislocation density of the matrix material side is of artificial 
nature due to an interpolation between the integration points. Due 
to the dislocation transport equation the total amount of dislocations 
is preserved inside the fiber-matrix composite. Consequently, due to 
the pile-up of dislocations at the fiber-matrix interfaces, the density of 
dislocations sinks in areas of homogeneous plastic slip inside the matrix 


material. 
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Figure 7.5: Illustration of the simple shear of a fiber-matrix composite material with 
periodic boundary conditions (top) and FE results shown in the zz-plane (bottom). The 
left-hand side shows the final micromorphic accumulated plastic slip jac. The right-hand 
side shows the corresponding coarse-grained accumulated density of dislocations pac. 
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7.7 Summary and conclusion 


An approach is presented which allows to model the transport of 
dislocations on a coarse-grained scale within a numerically efficient 
framework based on accumulated field variables. With the approach 
of Eq. (7.22) the in-plane gradient of the accumulated plastic slip is 
prorated on the individual slip systems which preserves the slip kine- 
matics of in-plane dislocation gliding. Additionally to the accumulated 
plastic slip the coarse-grained accumulated density of dislocations 
is introduced. With an energetic GB the dislocation pile-up at GBs 
(microhard) and the dislocation annihilation after activation of GB slip 
(microcontrolled) is modeled. This is consistent to the experimental 
work of Sun et al. (1998), where the buildup of dislocations near the 
GB and a vanishing peak at higher strain levels is observed. The FE 
implementation is discussed as well as the distributions of the plastic 
slip and dislocation density over a periodic laminate structure and a 
fiber-matrix composite material. The presented approach provides a 
promising numerically efficient framework for further investigations, 
for example investigations on the curvature induced line length pro- 
duction, which is neglected in this work. 

Regarding the continuity of plastic slip across GBs the following has 
to be considered: For the special case of equally oriented grains, as 
considered in this chapter, the accumulated plastic slip is continuous 
across GBs. As discussed in Chapter 5 for misoriented grains, however, 
a discontinuous solution field is required. One possible approach for 
misoriented grains is the enrichment of trial functions with a discontin- 
uous enrichment function as presented in Chapter 6. 
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Summary and outlook 


This thesis addresses several open questions of dislocation-GB inter- 
actions within gradient crystal plasticity frameworks. The behavior 
polycrystalline microstructures under consideration of dislocation pile- 
ups at GBs was realized by a model approach based on an extended 
energy balance according to Prahs and Bóhlke (2020a;b). An equiva- 
lence to the model approach based on an extended principle of virtual 
power according to Gurtin (2002) was obtained under the assumption 
of isothermal processes, a homogeneous distribution of the temperature 
and a temperature-independent Helmholtz free energy for the bulk ma- 
terial and the GBs. Under these assumptions, however, it was shown 
that both the bulk material as well as the GBs are dissipation free. 


The gradient crystal plasticity framework was applied to a three-phase 
periodic laminate microstructure in order to investigate dislocation 
transmission through GBs. The analytical solution of a boundary value 
problem was investigated for the special case of single slip. A GB 
energy according to Gurtin (2008) based on the geometric dislocation 
tensor was used in order to account for a GB resistance based on the 
grain misorientation angle For the limit case of a vanishing grain 
misorientation the single crystal consistency check is failed with the 
commonly used GB flow condition approach based on the assumption 
of independent plastic slips. In order to account for this shortcoming a 
new model approach for the GB flow condition was presented in this 
work. An additional flow rule is introduced in which the discontinuity 
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of the plastic slip gradient at the GB is correlated to a grain misorienta- 
tion function. 

Applying the new approach to the limit case of a vanishing grain mis- 
orientation results in two different solutions. While one solution is 
an artificial, non-physical solution identical to the solution obtained 
by the common approach, the second solution showed a single-crystal 
response with a continuous and symmetric parabolic distribution of the 
plastic slip even under variation of the grain size. The effect of the grain 
size on the mechanical response was investigated. A good agreement 
is found with a simple inverse or logarithmic approach, which is sug- 
gested in Li et al. (2016) based on underlying dislocation theories. 

With the introduced GB flow condition approach at first, a microhard 
state results in a pile-up of dislocations at the GB until a critical state 
is reached. The transition to a microcontrolled GB state where disloca- 
tion transmission across the GB occurs strongly depends on the misori- 
entation angle between adjacent grains. The presented approach pro- 
vides a promising framework the modeling of polycrystalline materials 
within a finite element implementation, where various misorientation 
angles and grain sizes are involved. For each slip system, however, one 
gradient-stress associated with the in-plane gradient of plastic slip has 
to be introduced. This results in a high computation cost, which is 
increased even more under consideration of the discontinuous distri- 
butions of the plastic slips which result from an orientation dependent 
dislocation transmission across GB. 

The approach of Erdle and Bóhlke (2017) was presented where one 
gradient-stress associated with the accumulated plastic slip is used as 
an overall measure of long range dislocation interactions within a finite 
strain setting. This provides a promising numerically efficient frame- 
work for large scale three-dimensional FE simulations. Discontinuous 
trial function were introduced in order to evaluate a GB flow rule in due 
consideration of a discontinuous distribution of the accumulated plastic 
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slip across GBs. The FE implementation of a periodic bicrystal was used 

to investigate the influence of energy parameters and size effects. 

An extension to the accumulated plastic slip framework was presented 

which allows to model the transport of dislocations on a coarse-grained 

scale. The in-plane gradient of the accumulated plastic slip is prorated 

on the individual slip systems which preserves the slip kinematics of 

in-plane dislocation gliding. A dislocation transport equation based on 

the coarse-grained accumulated density of dislocations was introduced. 

With the presented framework the two stages of GB behavior which are 

observed in Sun et al. (1998), i.e., buildup of dislocations near the GB 

and a vanishing peak at higher strain levels, was modeled based on 

energetic GBs. The dislocation pile-up at GBs (microhard) and the dis- 

location annihilation after activation of GB slip (microcontrolled) was 

observed for laminate microstructures and fiber-matrix composites. 

Even though multiple relevant results regarding dislocation-GB inter- 

actions have been obtained within this work, further investigations are 

necessary. In detail, these are: 

* Investigation of dislocation behavior at GBs under cyclic loading 

* Modeling of temperature dependent dislocation and GB effects 

e Comparison of gradient crystal plasticity framework based on indi- 
vidual slip gradients with the framework based on the gradient of 
the accumulated plastic slip 

* Modeling of an orientation-dependent GB flow rule in the accumu- 
lated plastic slip framework 

* Investigation of curvature induced line length production in the dis- 
location transport equation 

* Direct validation of dislocation-GB interactions with experimental 
investigations and discrete dislocation dynamics simulation 
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Appendix 


A.1 Slip system convention of face-centered 
cubic unit cell 


The slip system convention in reference configuration, which is used for 
the implementation of the FCC materials, is listed in Table A.1. The slip 
systems are labeled according to the Schmid-Boas slip system conven- 
tion, cf. Schmid and Boas (1935). 


Table A.1: Slip system convention for FCC in Schmid-Boas notation. 
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A Appendix 


A.2 Details on the global algorithmic tangent 


The total derivatives of Eqs. (6.47), (6.49) and (6.51) for the nodal DOF 
(à, 3, ö} read in matrix-vector notation 


T a an — O8 | fag 

ov x g 
a^|-|o$ 5 Flle! (A.1) 
dr? Et d AM 


The entries of the global algorithmic tangent are computed by 
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